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1  13  ABSTRACT 

This  is  Volume  I  of  a  finoi  report  presented  in  three  volumes:  Voi  I  -  Analysis  Methods;  Vol  !l  - 
Fabrication,  inspection  and  Testing;  Volume  III  -  Fatigue  Analysis  and  Failure  Mode  Studies. 

This  volume  describes  the  analytical  methods  that  were  generated  cr  used  as  part  of  th is  program. 
The  main  body  of  this  volume  is  divided  into  three  chapters  as  follows:  Closed  Foim  Analysis 
Methods,  Finite  Element  Anoiysis,  Photoelastic  Stress  Analysis.  Emphasis  wes  placed  on  rhe 
development  of  closed  form  analysis  procedures  for  bonded  joints  in  laminated  composites.  A 
comprehensive  lineai  analysis  method  and  associotea  computer  program  (KONJO  I)  has  been 
developed.  Numerical  resuirs  obteined  with  tMs  progrem  are  compared  wiln  finite  element 
analyses,  strain  gage  data  and  phatoeiostic  results  A  “Plastic  Zone”  approach  jvcs  used  to 
extend  BONJO  I  *o  include  joints  with  ideally  elastic. -plastic  adhesive  stress-strain  behavior. 

The  theoretical  development  of  a  rigorous  non-linear  analysis  procedure  fci  bonded  joints  has 
been  presented,  however,  this  method  was  not  carried  beyond  the  exploratory  stage. 

Finite  element  analyses  used  to  evoiuote  the  step  lap  bonded  join’s  and  bo'ted  joints  -'s  presented 
and  discussed.  Photoelastic  stress  analysis  procedures  used  in  the  program  ore  described  in  the 
final  chapter  together  with  the  results  obtained. 

Oetsiis  Of  iHustratiwi*  »n 
tt,*4  document  m ay  w  beU%f 

tiudi*}  ft?  IShTOtldO 


DD  .".r.,1473 


UNCLASSIFIED 

Stvytjh  CU4*.ncats**o 


UNCLASSIFIED 

Security  Classification 


key  words 


bonded  joints 
mechanical  joints 
joint  analysis 
boon  composite  materials 
fatigue  testing 
failure  modes 
fotigue  endurance 
fatigue  analysis 
photoelastic  stress  analysis 
mateiiaf  properties 
joint  fabrication 
non -destructive  inspection 


ROLE  I  #T  J  ROLE  1  WT  1  POL  b  J  W 


UNCLASSIFIED 

Vcuritv  Cl*»r'f:e»t-on 


DEVELOPMENT  OF  AN  UNDERSTANDING  OF  THE  FATIGUE  PHENOMENA 
OF  BONDED  AND  BOLTED  JOINTS 
IN  ADVANCED  FILAMENTARY  COMPOSITE  MATERIALS 

VOLUME  L  ANALYSIS  METHODS 


J.  N.  Dickson,  T.  M.  Hsu,  and  J.  M.  McKinney 
Lockheed -Georgia  Company 
Marietta,  Georgia 


Approved  fcr  public  release;  distribution  unlimited 


FOREWORD 


This  report  summarizes  the  work  accomplished  under  Contract  F33615-70-C-1302, 
"Development  of  an  Understanding  of  the  Fatigue  Phenomena  of  Bonded  and 
Bolted  Joints  in  Advanced  Filamentary  Composite  Materials",  and  was  prepared 
by  the  Lockheed-Georgia  Company,  a  Division  of  Lockheed  Aircraft  Corporation. 
The  work  reported  herein  was  sponsored  by  the  Advanced  Composite  Branch,  Air 
Force  Flight  Dynamics  Laboratory,  Air  Force  Systems  Command,  Wright-Patterson 
Air  Force  Base,  Ohio  45433.  Mr.  Rodman  Joblove,  FBC,  was  the  Air  Force 
Project  Engineer  and  Mr.  A.  C.  Fehrle  was  the  Lockheed-Georgia  Program 
Manager. 

The  authors  of  Volume  I  are  Mr.  J.  N.  Dickson,  Dr.  T.  M.  Hsu,  and 
Dr.  J.  M.  McKinney.  Mr.  Dickson  and  Dr.  Hsu  are  jointly  responsible  for 
the  development  of  the  ..losed  form  analysis  procedures  and  the  finite  element 
analyses  performed  during  the  program.  Dr.  McKinney  was  responsible  for  the 
photoelastic  analyses. 

This  technical  report  has  been  reviewed  and  is  approved. 

For  internal  control  purposes,  this  report  has  been  assigned  Lockheed-Georjia 
Company  Report  Number  ER-1 1319. 


Philip  A.  Parmley  £ 

Chief,  Advanced  Composite  Branch 
Structures  Division 

Air  Force  Flight  Dynamics  Laboratory 


ABSTRACT 


This  volume  describes  the  analytical  methods  thot  were  generated  or  used  os  port  of  this 
program.  The  main  body  of  this  volume  is  divided  into  three  chapters,  as  follows: 

o  C:osed  Form  Analysis  Methods 

o  Finite  Element  Analyses 

o  Photoelastic  Stress  Analysis 

Emphasis  wcs  placed  or<  the  development  of  closed  form  analysis  procedures  for  bonded 
joints  in  laminated  composites.  A  comprehensive  linear  onclysis  method  and  associated 
computer  program  (BONJO  I)  has  been  developed.  Numerical  results  obtained  with  this 
program  are  compared  with  finite  element  analyses,  strain  gage  data,  and  photoelastic 
results.  A  "Plastic  Zone"  approach  was  used  to  extend  BONJO  I  to  include  joints  with 
ideally  elastic -plastic  cdbesive  stress-strain  behavior.  The  theoretical  development  of  a 
rigorous  non-iinsar  analysis  procedure  for  bonded  joints  has  beer  presented.  However, 
this  method  was  not  carried  beyond  the  exploratory  stage. 

Finite  element  analyses  used  to  evalucre  the  step  lap  bonded  joints  arid  bolted  joints  are 
presented  and  discussed.  Photoelostic  stress  analysis  procedures  used  in  the  program  are 
described  in  the  final  chapter  together  with  the  results  obtained. 
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I.  INTRODUCTION 


This  program  was  undertaken  to  develop  an  understanding  of  the  fatigue  phenomena  of 
structural  joints  in  advanced  filamentary  composite  materials  and  to  develop  analytical 
and  testing  methods  to  support  proper  fatigue  design  of  advanced  composite  structural 
joints.  The  program  included  the  evaluation  of  both  bonded  and  bolted  joints.  Primary 
emphasis  was  placed  on  joints  in  boron-epoxy.  However,  a  limited  evaluation  of  bonded 
joints  in  graphite-epoxy  and  glass-epoxy  was  included.  Although  *\e  sizes  of  the  joints 
for  this  investigation  were  small  (one  to  ten  inches  in  width)  all  configurations  evaluated 
are  representative  of  typical  structural  joints  currently  utilized  in  advanced  filar ientary 
composite  structures. 

The  program  consisted  of  three  major  areas  of  investigation: 
o  Analysis  Methods 
o  Fabrication,  Inspection,  and  Testing 
o  Fctigue  Analysis  and  Failure  Mode  Studies 

Analytical  methods  for  determining  joint  stresses  were  divided  into  two  major  tasks: 

(1)  analysis  of  bonded  joints  and  (2)  analysis  of  bolted  joints.  Primary  emphasis  was 
placed  on  the  development  of  a  closed  form  elastic  analysis  procedure  for  bonded  joints. 
This  analysis  was  used  to  evaluate  a  number  of  joint  variables.  A  "plastic  zone"  approach 
was  used  to  extend  the  closed  form  analysis  procedure  to  include  joints  with  inelastic 
adhesive  stress-strain  behavior.  The  results  of  the  elastic  closed  form  solution  were 
verified  with  finite  element  analyses,  photoelastic  analysis  and  strain  gage  data.  Finite 
element  analyses  were  used  to  evaluate  the  step  lap  bonded  joints  and  bolted  joints. 

The  experimental  program  consisted  of  fabrication,  inspection  and  testing  of  a  large 
quantity  of  joint  specimens.  Fabrication  and  inspection  methods  were  established  which 
resulted  in  specimens  being  fabricated  to  close  tolerances  and  of  uniformly  high  quality. 
This  provided  specimens  that  would  consistently  develop  stresses  that  were  predicted  by  the 
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analytical  methods.  Developing  testing  techniques  and  actual  specimen  testing  was  a 
major  portion  of  the  program.  Establishing  proper  specimen  support  was  essential  to 
obtaining  repeatable  joint  strengths  within  a  specimen  configuration.  Equally  important 
was  determining  the  proper  cyclic  rate  for  the  different  stress  ratios  and  specimen  configu¬ 
rations  to  preclude  specimen  heating  and  erratic  fatigue  lives. 


Evaluation  of  the  experimental  results  was  divided  into  two  seprrate  but  related  tasks. 
These  tasks  were  failure  mode  studies  and  fatigue  analysis.  The  failure  mode  studies 
mentioned  were  photomicrographic  analyses  of  the  failure  surfaces.  This  failure  mode 
analysis  does  not  replace  but  augments  the  gross  failure  modes  generally  defined  within  the 
experimental  phases  of  a  program.  The  photomicrographic  analysis  conducted  within  this 
program  established  failure  modes  related  to  specific  joint  designs,  joint  loading,  and 
fatigue  history.  The  fatigue  analysis  established  relationships  between  specimen  configu¬ 
rations,  joint  variables,  material  combinations,  loading  conditions  and  stress  ratio  effects 
for  constant  amplitude  loading.  The  relationship  between  constant  amplitude  fatigue  and 
spectrum  fatigue  (block  and  realistic)  was  also  evaluated  for  specific  joint  configurations. 


This  report  is  divided  into  three  separate  volumes  each  containing  the  developments 
accomplished  within  a  major  area  of  investigation.  Each  volume  is  a  self-contained 
document,  complementing  the  other  two  volumes  but  not  dependent  upon  them  for 
coherence  or  continuity.  The  titles  of  the  three  volumes  are: 

Volume  I  -  Ana  lysis  Methods 

Volume  II  -Fabrication,  Inspection  and  Testing 

Volume  III  -  Fatigue  Analysis  and  Failure  Mode  Studies 


Volume  I  is  divided  into  three  chapters,  with  primary  emphasis  being  placed  on  the  first 
chapter.  This  chapter,  entitled  "Closed  Form  Analysis  Methods,"  contains  sections  on 
earlier  methods  and  their  short-comings;  the  development  of  the  analysis  procedure 
BONJO  I;  the  plastic  zone  extension  of  BONJO  I;  and  a  section  on  numerical  results. 
The  second  chapter  is  devoted  to  finite  element  analyses  of  the  step  lap  bonded  joints  and 
mechanical  joint  specimens  while  the  third  chapter  deals  with  photoelastic  analysis 
methods. 
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I  •  I 

II.  GENERAL  DISCUSSION 


1.0  INTRODUCTORY  REMARKS 


The  increasing  utilization  of  advanced  filamentary  composite^  in  the  design  of  aircraft  and 
space  vehicles  has  necessitated  the  development  of  more  sophisticated  methods  of  analysis. 
Prior  to  the  introduction  of  these  new  mqteriaisf,  the  use  of  adhesive  bonding  was  limited 
almost  entirely  to  secondary  structural  applications.  The  design  of  bonded  joints  was  based 
largely  on  experimental  data.  Theoretical  methods  such  as  those  by  Goland  and 

i 

Reissner  (1),  Volkerson  (2),  and  Szepe  (3)  were  used  to  obtain  a  more  detailed  knowledge 
of  the  distribution  of  adhesive  shear  and  normal  stresses  in  bonded  lap  joints  and  were 
generally  considered  as  adequate  for  that  purpose. 


The  determination  of  stresses  in  bonded  joints  with  laminated  composite  adherends,  . 
however,  is  much  more  complex  than  In  those  with,  isotropic  adherends.  This  increased 
complexity  can  be  only  partly.atfributed  to  the  fact  that  composites  are  anisotropic  and 
heterogeneous.  A  large  part  is  caused  by  the  necessity  to  account  for  the  effects  of  inter¬ 
laminar  shear  and  normal  stresses  in  the  analysis  of  composite  joints.  Numerical  results 
that  demonstrate  the  importance  of  the  above  effects  are  presented  in  tnis  report. 

s 

Another  factor  that  may  be  important  in  the  analysis  of  composite  joints;is  the  presence  of 
residual  thermal  stresses.  These  stresses  ar'e  caused  by  bonding  or  curing  at  elevated 
temperatures  and  subsequent  cooling  to  operating  temperatures.  Although  thermal  stresses 
are  introduced  in  all  multi-directional  laminates,  they  become  especially  significant  when 
bonding  highly  dissimilar  materials  such  as  boron  and  aluminum  or  graphite  and  aluminum. 

It  is  within  the  capability  of  most  large  finite  element  programs  in  existence  today  to 
perform  an  analysis  of  virtually  any  degree  of  complexity;  but  such  analyses  are  cumber¬ 
some  and  costly.  In  addition,  the  preparation  of  input  data  and  interpretation  of  results  is 
extremely  time  consuming.  It  is  therefore  preferable  to  use  direct  or  closed  form  methods 
wherever  possible. 
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Since  the  majority  of  specimens  to  be  analyzed  as  part  of  this  program  were  simple  bonded 
lap  joints  of  uniform  geometry,  it  was  decided  to  develop  a  closed  form  procedure  for 
composite  joints.  Although  it  became  necessary  to  develop  a  more  comprehensive  pro¬ 
cedure  than  initially  anticipated,  for  the  reasons  mentioned  previously,  the  method  turned 
'out  to  be  very  efficient  and  offers  tremendous  cos.  and  time  savings  as  compared  to  the 
finite  element  method.  Finite  element  analyses,  however,  were  performed  for  some  of  the 
step  lap  joint  and  mechanical  joint  specimens  and  for  the  verification  of  the  closed  form 
procedure.  ^ 


2.0  NOMENCLATURE 


Symbols  are  defined  in  the  text  when  they  first  appear.  A  list  of  the  most  important  ones 

i 

is  given  below  for  convenience. 

A,  B,  D  Laminate  Stiffness  Matrices 

E  .  Young's  Modulus 

G  Shear  Modulus 

L  Overlap  Length 


Stress  Couples 
Stress  Resultants 
Transverse  Shear  Forces 

Distance  from  Reference  Surface  to  Free  Surface  of  Laminate 
Thickness 


u,  v,  w  Displacements  in  x,  y,  and  z  directions 

x,  y,  z  Cartesian  Coordinates 
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a  ,  a 

X  y 


£  /  ^  /  c 

x  y  xy 


£xp'  €yp'  €xyp 


Thermal  Expansion  Coefficients 


In  Plane  Strains 


In  Plane  PUjstic  Strains 


Normal  Strain 


zp 


e  ,  e 
xz  yi 


Normal  Plastic  Strain 
Transverse  Shear  Strains 
Adhesive  Effective  Stress 


a  ,  c  ,  a 
x  y  xy 


In  Plane  Stresses 


Normal  Stress 


:  ,  a 

xz  yz 


v 


Transverse  Shear  Stresses 


Poisson's  Ratio 


u;  L,  a 
Superscript 


Upper,  Lower,  Adhesive 


th 


Reference  Surface,  i<  Lamina,  Temperature 
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ill.  CLOSED  FORM  ANALYSIS  METHODS 


1.0  GENERAL 

For  the  purpose  of  this  report,  theoretical  methods  other  than  finite  element  or  finite 
difference  procedures  will  be  defined  as  c  losed  form  methods  of  analysis.  This  will 
include  ite-ative  as  well  as  direct  solutions  based  on  the  theory  of  elasticity.  The  assump¬ 
tions  of  small  deformations  and  uniform  geometry  are  made  throughout  this  section. 
Additional  assumptions  will  be  discussed  as  they  are  applied.  Closed  form  solutions  are 
developed  for  the  two  joint  configurations  shown  in  Figure  1 .  Because  of  assumed  sym¬ 
metry,  only  one  half  of  the  single  lap  joint  and  one  quarter  of  the  double  lap  joint  need  to 
be  considered  as  shown  by  the  heavy  dashed  line.  The  length  of  the  joint  is  defined  in 
Figure  1  as  being  in  the  plane  of  the  paper  and  in  the  direction  of  the  chosen  x-csxis.  The 
width  of  the  joint  is  the  dimension  perpendicular  to  the  plane  of  the  paper  (parallel  to  the 
y-axis} .  The  region  OsxS  L  will  be  of  primary  interest  and  the  governing  differentia! 
equations  to  be  developed  will  apply  to  this  region  only. 


2.0  PRELIMINARY  INVESTIGATIONS 

After  a  thorough  evaluation  o'  available  theoretical  methods,  it  was  decided  early  in  the 
program  to  develop  a  procedure  for  the  anu  lysis  of  simple  laminated  composite  lap  joints 
based  cn  the  Goland-Reissner  (I)  differential  equation  approach.  This  procedure  required 
the  removal  of  the  Imposed  symmetry  restrictions  in  the  Goland-Reissner  analysis  and  the 
substitution  of  gross  laminate  properties  in  p.lace  of  the  isotropic  elastic  moduli.  The 
K ire hhoff assumptions,  however,  were  assumed  to  remain  valid  in  this  initial  procedure. 

For  the  case  of  single  lap  joints  with  relatively  flexible  bond  layers  and  isotropic 

adherends,  Goia.nd  and  Reissner  obtained  the  following  uncoupled  differential  equations 

for  the  shear  stress,  c°  ,  and  the  normal  stress,  a°,  in  the  adhesive: 
xz  z 
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DETAIL  A 


FIGURE  1  LAP  JOINT  CONFIGURATIONS 
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It  might  be  remarked  that  there  is  an  inconsistency  in  equations  (i) ,  The  first  equation  is 
based  on  plane  stress  considerations  while  the  second  one  assumes  a  condition  of  plane 
strain. 

The  above  equations  are  uncoupled  as  a  result  of  the  requirements  that  the  thickness  and 
material  properties  of  the  two  adherends  must  be  identical .  If  these  restrictions  are 
removed  one  may,  by  using  essentially  the  same  procedure,  obtain  the  two  coupled  dif¬ 
ferential  equations: 


d3a°  4G 
xz  a 

/i-,2 

1  U  + 

i  2' 

1 

dx3 

t 

a 

\E  t 
\  u  u 

ElA> 

4  o 

da  12E  j 
z  +  °l 

M+ 

1 

dx4 

'a  1 

W 

6Ga/' 


E  t2 
u  u 


CT°  =0 
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la°  +  T 


u  2 

1  -V 
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1 

Ida 

xz 

c  2 

c  2  , 

f  dx 

.  E  t 
^  u  u 

W) 

These  equations  thus  represent  the  governing  differential  equations  for  a  single  lap  joint 

with  unequal  isotropic  adherends  in  a  state  of  plane  strain  (in  the  x-z  plane).  E^  and  E^ 

are  the  elastic  moduli,  v  and  v.  are  the  Poisson’s  ratios,  and  t  and  t.  are  the  thicknesses 

u  L  u  L 

of  the  upper  and  lower  adherends,  respectively.  To  obtain  the  relations  for  a  plane  stress 
condition,  the  Poisson's  ratios  are  dropped  from  the  above  equations. 


The  norma!  stress,  a  ,  may  be  found  from  the  first  of  equations  (2)  and  then  substituted 
into  the  second  one  to  yield  the  following  sixth  order  differential  equation  ir*  (da^^/dx): 

/  A®  j4  j2  \  da° 

(-^T+C,  ~+c9^T  +  cJ-^  =  °  (: 

\dx6  1  dx4  2dx  3/  dx 
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In  the  above  represents  the  total  thickness  of  the  plates  to  be  joined  and  t  is  the  thick¬ 
ness  of  each  splice  plate. 


% 

I 


The  differential  equation  (3)  applies  to  joints  with  laminated  composite  adherends;  but  the 
coefficients  C.  will  now  involve  the  gross  laminate  stiffness  properties.  The  stiffnesses 
may  be  obtained  from  laminated  plate  theory*  without  difficulty.  The  constitutive  equa¬ 
tions  for  a  thin  laminated  plate  are  generally  written  in  the  form: 

{N}  =  [A](e]  -  [B]{k] 

{M}  =[B]{e}  -  [D]{k] 

When  the  laminate  is  in  a  state  of  plane  strain,  one  has,  therefore: 


(7) 


N  =  A,  ,e 
x  It  x  tlx 


M  =  B,  ,e  -D..x 

X  tlx  1  1  X 


(3) 


where  e°  is  the  strain  at  the  reference  surface,  which  will  for  convenience  be  taken  at  the 
x 

interface  of  the  laminate  and  the  bond  layer.  The  gross  laminate  stiffnesses  are  given  by: 

N 

An'±Icn>Vhk-i> 

k=l 

BU  *42  C'n  )(hk  "  hk-i) 

k=l 


(9) 


N 


n  .  1v  r‘(k)/L 3  ,3  . 

Dn  3  Z  cn  *hk " hk-i* 

k=l 


l 


Dn“Dli  -  B1VA11 


*Note  that  at  this  point  the  Kirchhoff assumptions  are  still  assumed  to  be  valid. 
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W 
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H 
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In  the  above  expressions,  the  positive  sign  refers  to  the  upper  laminate  and  the  negative 

th 

sign  to  the  lower  laminate;  h^  is  the  distance  from  the  upper  surface  of  the  k  layer 
(lower  surface  for  lower  laminate)  to  the  bond  layer  interface.  For  a  single  iap  joint  with 
composite  adherends,  one  now  obtains: 


C  =-S|_L+-L 

2  t  Uu  rL 

a'°ll  Dll, 


C3=ClV 


E  G 
a  a 


B 


11 


B 


11 


Uu  DU  AL  D1" 

VAnun  aiiuil 


while  for  a  double  lap  joint  one  finds 

G 


D 


c,= 


n 


t 


K,  M 


E  , 

M 

— 

t 

a  ' 

(°u  / 

C3  C1C2  + 


E  G 
a  a  i 


B 
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2  1  .U  prU 

a  VAn  D1 1 


(4c) 


(4d) 


£ 


When  the  laminates  are  in  a  state  of  plane  stress  (in  x-z  plane),  i.e. 


oW-o«-„«-0 

y  xy  yz 

'(h)  (k) 

the  gross  laminate  stiffness  properties  are  obtained  by  replacing  C ^  '  by  1/5^  j  in  equa¬ 
tions  (9).  and  Sj  ^  are  defined  in  reference  4  in  terms  of  the  elastic  constants  and 
lamina  orientations. 
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After  the  coefficients,  C.,  are  determined,  the  roots,  X.,  can  be  calculated  from  equation 
(6).  Seven  boundary  conditions  are  now  required  to  obtain  a  solution  and  define  the  state 
of  stress  in  the  joint.  None  of  the  conditions  that  were  used,  however,  pertained  to  the 
adhesive  shear  stress  and  as  a  result  the  peak  shear  stresses  occurred  at  the  ends  of  the 
joint,  as  was  the  case  in  the  Goland-Reissner  analysis.  This  cannot  be  correct,  of  course, 
and  it  will  be  shown  later  that  in  reality  the  peak  stresses  occur  slightly  inward  from  the 
edges. 

The  fact  that  the  location  of  the  peak  shear  stresses  in  the  adhesive  could  not  be 
accurately  predicted  was  not  considered  a  major  argument  against  the  use  of  the  method 
described  here,  but  studies  were  made  which  indicated  that  the  magnitude  of  these  stresses 
was  considerably  lower  than  calculated  with  this  analysis.  It  was  conjectured,  and  later 
substantiated,  that  because  of  the  relatively  low  transverse  shear  stiffness  and  normal 
stiffness  (through  the  thickness  of  the  laminate)  of  most  fibrous  composites,  adherence  to 
the  Kirchhoff  assumptions  could  lead  to  unacceptable  errors. 

The  above  method  of  analysis  was  therefore  discontinued  and  a  more  rigorous  procedure 
developed  in  which  the  Kirchhoff  assumptions  were  discarded.  The  theoretical  development 
of  this  procedure  is  presented  in  the  following  sections. 

3.0  GENERAL  EQUATIONS  FOR  LAMINATED  COMPOSITE  JOINTS 

3.1  Lamina  Relations 

The  laminate  is  assumed  to  consist  of  orthotropic  layers  which  may  have  arbitrary  thickness 
and  material  properties.  A  Cartesian  (x,  y,  z)  coordinate  system  is  used  as  the  plate 
reference  system  in  which  the  z-axis  is  normal  to  the  plane  of  the  plate.  A  small  splice 
cut  from  the  laminated  plate  is  shown  in  Figure  2.  The  material  symmetry  axes  (1 ,  2)  of 
each  lamina  lie  in  a  plane  parallel  to  the  x-y  plane  but  may  be  rotated  about  the  z-axis 
through  an  arbitrary  angle  8.  The  thermal  expansion  coefficients  are  also  assumed  to  be 
orthotropic  and  hence: 
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«'V -'  “  *  W  ■WJWWI^ 


AS'fn,  a  p*j ^pip-i!j«ieg,*r5w5a  i 


a4  =  a5  =  a6  =  0 


th 


The  linear  strain-displacement  relations  for  the  k  lamina  may  be  written: 


.00 

=  u(k) 

e(k)  = 

= VW  + 

w(k) 

X 

u  / 

X 

yz 

V'z 

'y 

.00 

-;k> 

e<k>  = 

(k) 

w. 

y 

y 

xz 

z 

'x 

(k) 

-W« 

e(k)  = 

=  0W  + 

v(k) 

z 

'z 

xy 

'y 

V'x 

(10) 


where  a  comma  foiiowed  by  a  symbol  means  differentiation  with  respect  to  that  symbol. 

th 

The  constitutive  relations  for  the  k  lamina  with  respect  to  its  own  principal  (1 ,  2,  3) 
axes  may  be  expressed  in  matrix  form  by 


r0"i(k)  =  m(k)  fc-  nil  (k) 

0  1,2,3  'C^l ,2,3  ‘aT^  1,2,3 


(ID 


T  is  the  temperature  change  from  some  initial  reference  state  (bonding  or  curing  tempera¬ 
ture,  for  instance) .  The  non-zero  coefficients  of  the  stiffness  matrix  [C]  are  given  in 
reference  4  in  terms  of  the  elastic  constants  of  the  material .  The  stress-strain  equations 
with  respect  to  the  laminate  reference  (x,y,z)  axes  are: 

(k) 


' _ 

a 

X 

(k) 

Q11  Q12  Q13  0  0  Q16 

(k) 

'  > 

e 

X 

a 

y 

Q22  Q23  0  0  Q26 

e 

y 

a 

z 

► 

Q33  0  0  Q36 

e 

z 

a 

o 

i 

O 

a 

z 

yz 

44  45 

yz 

a 

0 

1 

t 

o 

G 

xz 

55 

XZ 

a 

Qz, 

e 

xyj 

66 

xyj 

(12) 


Since  the  [Q]  matrix  is  symmetric,  only  the  upper  half  of  the  matrix  is  shown.  The  actual 
stress  in  the  lamina  is  obtained  from  the  relation 


{a}(k)  ={a}(k)  -  £aT}(k>  =  [Q](k){el(k>  -  {cT}(k) 

x,y,z  Jx,y,z  Jx,y,z  1  'x,y,z  1  jx,y, 

where  the  equivalent  thermal  stresses  aie  defined  cs  follows 


(13) 
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[(CUal  +C12tt2  +  C13a3)m2+  <C12al  +  C2ft  +  C23a3)n2;1T 


[(C)  lal  +  C12a2  +  C13a3)n2  +  (C1^1  +  C2^2  +  C23a3>m2lT 


(Cl3al  +  C23a2  +  C33a3)T 


T  _  T  - 
a  =  a  =0 
yz  xz 


°xy  =  l‘Cllal  +  Cl2(a2  -al>  -C2*+(C13  •C23)a3]mnT 


and  m  =  cos  8,  n  =  sin  0.  The  Q..  are  given  by: 


4  2  2  4  2  2 

Qjj  =  m  C^  +  2m  n  +  n  ^22  +  ^fn  n 


Q^2  =  m^n^C^  +  (m^  +  + 


Q13  =  m  C13  +  n  C23 


Q16  =  m^nC^  ^  -  mn(m^  -  n^JC^  ~  mn*^-22  ~  2mn(m^  -  n^)C^ 


4  2  2  4  2  2 

^22  =  n  ^*11  +  ^m  n  ^12+m  ^22  +  n  ^“66 


Q23  nC13  +  mC23 


^26  ""  mn^*]  ]  *  mn(m^  "  n^*]2  ~  m^n^*22  +  ~  n^66 


Q33  C33 


Q36  mn^C13~C23^ 


Q44  mC44+nC55 
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''7*  /srp&pyem ■sEf  £?>**»  «v~  ,-,* 


>■»*  -  •+=»**►  -* 


Q45='"’nC44+"“’C55 


Q55  =  n2C44+n,2c55 


~  _  2  2r  ,  2  2r  ,  2  2r  ,  ,  2  2.2  r 

®66  _m  n  ^*11  -  2m  n  C12+m  n  C22  +  m  ~  n  ^*66 


3.2  Effect  of  Transverse  Shear 

In  lap  joints  made  of  conventional  metallic  materials  the  effect  of  transverse  shear 
deformation  on  the  stress  distributions  in  the  joint  is  generally  found  to  be  negligible 
and  therefore  an  analysis  based  on  the  classical  Kirchhoff  assumptions  is  adequate  in  most 
cases.  However,  in  laminated  plates  having  relatively  low  transverse  shear  moduli,  these 
assumptions  are  no  longer  permissible.  The  effect  of  transverse  shear  deformation  on  the 
bending  of  symmetrically  laminated  plates  has  been  investigated  by  Whitney  (5) .  His 
work  has  been  extended  for  use  in  the  present  analysis. 


The  mid-surface  of  the  adhesive  layer  has  been  taken  as  the  reference  surface  for  both 

laminated  plates.  The  adhesive  shear  stresses  at  the  reference  surface  are  denoted  by 

a°  and  a°  and  the  transverse  shear  stresses  in  the  k*  layer  are  assumed  in  the  form: 
xz  yz  ' 

a*l  =  “xz  ('  -  h^) J +  CQ55f(z)  +  b55z  +  -mV  +  CQ45f(z)  +  b45z  +  a45]tpy 
°yz  =  V  ('  -  i^)  +  tQ“f<z>  +  b45z  +  «4S*x  +  CQS«Z>  +  b44z  +  °44\ 


05) 


where  cp  and  cp  are  functions  of  x  and  y  only;  z  is  measured  from  the  reference  surface; 
x  y 

and  hence  will  be  negative  for  the  lower  laminate  (see  Figure  1).  It  must  be  noted  that  in 
this  case  h^  will  also  be  a  negative  quantity.  In  an  isotropic  material  the  transverse 
shear  stresses  vary  paraboliccliy  through  its  thickness.  Whitney's  assumption  of  a  para¬ 
bolic  variation  of  shear  stress  within  each  layer  therefore  seems  reasonable.  With  this 
assumption,  f(z)  is  taken  in  the  form: 


06) 


-7-,  s- ...... 


f(z)  =  -4-  (ihN  -  I2)  • 
hN 

which  is  symmetric  with  respect  to  the  middle  surface  of  the  laminated  plate  and  goes  to 

(k)  ' 

zero  at  the  outside  surfaces.  The  constants  a...  are  determined  from  the  requirement  that 

•I  • 

the  shear  stresses  are  continuous  at  the  interfaces  of  adjacent  layers,  or 

aJ^^aJ^  +  CQ^  -q!!*+1)]  f(hk)  i, j  =  4,5  k  =  1,2, ..  .(N-l)  (17) 


It  follows  from  equations  (15)  that: 


■.!!>.  o 


(k) 

and  hence  the  remaining  a).7  can  be  calculated  from  equation  (17).  Since  f(0)  =  f(h^)  =  0 
one  finds  for  the  N^1  lamina  !  i 


The  b..  follow  from  the  condition  that  the  shear  stresses  ore  zero  at  the  free  surface  . 
•I 

(z  =  h^)  and  therefore: 

(N) 

a.. 

b..  =  --jiJ —  (l; 

'i  hK1 

'•  N  i  i 


For  a  laminate  which  is  symmetrical  with  respect  to  its  middle  surface,  a..  7  and  therefore 

b..  will  be  zero.  1 

M 

I 

The  transverse  shear  forces  are  obtained  by  integrating  the  shear  stresses  over  the  thickness 


of  the  laminate 


±Vx=f°°zl,N+A5^x  +  A45cp, 


±V  =7 hK1  +  A  .Ccp-  +  A .  .ce 
y  2  yz  N  45^x  44vy 
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where: 


\  ■  * tyv.  +  I  tQ!^F\>  -  F\->» + 

k=l 

*/i  =  4/5 

and 

F(z)  =ffU)  dz  =  (3zhN  ■  2Z2) 

3hN 

i 

In  equations  (19),  the  positive  sign  refers  to  the  upper  laminate  or  splice  plate  while  the 
negative  sign  refers  to  the  lower  laminate  (negative  z-axis).  This  convention  will  be 

s  * 

maintained  throughout  the  cnaiysis.  Solving  for  the  shear  functions  cp  and  cp  from 

x  y 

equations  (19)- one  obtains: 

cd  =A_(±V  -ihMa°  )  +  A.C(±V  -ik.aC  ) 

x  55  \  x  2  N  xz/  45  \  y  2  d  yz/ 

cp  =  A  (±V  -ihK.a°  )+A  .(±V  -jrhMa°  \ 

y  45  \  x  2  N  xz /  44  \  y  2  N  yz  / 


fA-. 

A 

A 

A-  1 

1  44 

1 

45 

= 

44 

45 

l  45 

LO 

.  10 
l< 

-A45 

A55. 

Substituting  (2G)  into  (15/  yields  for  the  transverse  shear  stresses  in  the  k  layer: 

»_p(k)  o  ,p<k)  o  «<k)v  ^k),, 

°xz  “G1  axz+G2  V*  3  Vx±G.l  Vy 

,(k)=G«a°  +  G<kV  ±G<k>V  ±G<k»V 
yz  5  xz  6  yz  /  x  8  y 

(k) 

in  which  the  Gj  's  are  functions  of  z  only-  They  are  given  by: 

Gl<2)(k>=|-i^-5hNG3(='(ki 
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~  ,  s(k)  _  ] .  F  ,JL  k) 

G2  z  2  G4  2 

G3(z)(k)  =  VQ^ffzi  4  b,5z  v  o«]  i-  A4s[Q«fU>  4  b45z  ♦  0<;>1 

+  b45z  ♦  a«3  4  A45[Q<k5>f(2)  4  bMz  4  .«J 

S5(2)lk>==-|hN^2>W 

g6(z)(  *  -  1  -  -  2  hN  g8(z)1  ' 

S7U)(k)  =  A55[Q^f(z)  4  b45z  +  a®]  +  A45[Q«f(z)  +  b^z  4  0‘k)] 
58(z)(k)  -  A^QWffz)  4  kuz  4  .«!  4  A45[Q<k)f(z)  4  b43Z  4  .«] 


The  transverse  shear  strains  are  obtained  *rom  equations  (12): 


yz 

(k) 

Q44  Q45 

(k)  ' 

a 

yz 

(k) 

”S  S  T 

44  45 

xz 

°45  Q55. 

a 

xz 

c  c 

Lb45  "55J 

(k) 


yz 


xz 


(k) 


(22) 


ana  by  substitution  of  equations  (21)  one  has: 

(k)  _  rc(kMk)  .(kWk),  O  ,  r,(kWk)  (kW(k)  o 
exz'CS5SGl  S45G5  J°xz+C555G2  +  S45G6  ]ayz 

i  £S<k>5<k>  4  S»V  i  [S(kMk)  +  S<k»G(khv 

55  3  45  7  x  55  4  45  8  y 

e(k)  =  [s^k)G(,k)  +  S^G(k)]a°  +  [S^k)G(k)  + 

yz  45  1  44  5  xz  45  2  44  6  yz 

±  [S^G^  +  S^k)G(k)]V  ±  [sJJgJ*  +  S^G^IV 

4j  3  44  7  x  45  4  44  8  y 

The  in-plane  displacements  may  be  obtained  by  integrating  the  shear  strains  with  respect 
to  z.  From  equations  (10)  and  (23)  one  obtains,  therefore: 


(23) 
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(k)  o  f  ■  ,  _(k)  o  0(k)  o  r  (k). .  _  (k). . 

u  =  u  -  /  w,  dn  +  G,  a  +  GX  a  ±  GX  'V  *  O',  V 

J  'x.  1  xz  2  yz  3  x  4 


(k)  „  o  f 

v  —  v  -  : 


/  w,  dr>  *  G(k)a°  +  G^'a°  ±  G(-.k)V  ±  G(k)V 
J  y  5  xz  6  yz  /  x  8 


(24) 


where: 


G.(i)(k'=  /  [S*^5 .(n)®4  +  S45<3.+4W<kh*l  ■*  c| 


(« 


Gi+4(>),k)  -  /  ts^g-M00  +  +  .» 

°  i  =  1,2,3, 4 

Continuity  of  displacements  at  the  interfaces  of  adjacent  layers  requires  that: 


(25) 


G.(hk)(k)  =  G.(hk)(k+1)  (26) 

From  equations  (24)  it  can  be  seen  that  G.(0)^,  end  therefore  c^}\  is  zero.  The 

(k)  *  * 

remaining  constants  of  integration,  cj  ,  may  now  be  calculated  from  equation  (26). 

The  in-plane  strains  are  obtained  from  equations  (10)  by  differentiation  of  the  displace¬ 
ments  (24) .  Hence 


(k)  o 
e  -  e 


C  J-  a-  /~(k)  °  t  o(k)  °  r*(k)v/  /~(k)\/ 

-  /  w,  dtj  +  G,  a  +  GX  a  ±  GX  V  *  GX  V 

J  xx  l  xz,x  2  yz,x  3  x,x  4  y,) 

o 

-  f  w,  dri+G^a0  +  G*k)a°  ±  G^V  ±  G^V 

J  yy  5  xz,y  6  yz,y  7  x,y  8  y,) 


(27) 
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i 


H 


I 


>)  -  r° 

6  —  £ 

xy  xy 


-2  f  w,  dq  +  G(k)a°  +G(,k)a°  +  G^a° 

./  xy  5  xz,x  1  xz,y  6  yz, 

o 

i-  G(k)o 


(27) 

(cont'd) 


,  u°  ±  G^V  *  cikV  ±  G(k)V  ±  Gj.k)V 
2  yz,y  7  x,x  3  x,y  8  y,x  9  y,y 


Ak) 


The  functions  G.(z)1  '  thus  represent  the  effect  of  shear  deformation  in  the  laminate.  If  it 

is  assumed  that  the  normal  displacement,  w,  does  not  vary  through  the  thickness  of  the 

laminate,  the  integrals  :n  equations  (27)  ere  replaced  by  zw°  ;  zw°  ;  and  2zw°  ; 

xx  yy  ,xy 

respectively.  This  assumption,  however,  is  not  necessary  for  the  special  cases  considered 
in  this  report. 


3 . 3  Laminate  Stress  Resultants  and  Stress  Couples 


It  will  be  convenient  to  eliminate  the  normal  strain  from  the  lamina  stress -strain  relations 

th 

and  express  the  stresses  in  the  plane  of  the  k  lamina  in  terms  of  the  in -plane  strains  and 
(k) 

the  noimal  stress,  .  The  normal  strain  follows  from  the  third  of  equations  (12): 

(k)  _  1  frvk)  Q(k)  (k)  (k)  (k)  0(k)  (k)1 

ez  _^Jk)  L  z  -Q13ex  -Q23€y  'Q36ExyJ 


(28) 
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After  substitution  of  the  above  in  the  first,  second,  and  sixth  of  equations  (12),  the  in-plane 
stresses  may  be  written  in  the  form: 

(k)  f„.  l(k) 


a 

x 

(k) 

'Q'.i 

Q12 

*.6 

(k) 

e 

X 

(k) 

Q13 

a 

y 

►  = 

QI2 

Q22 

Q26 

£ 

y 

•  +' 

Q23 

o  . 

xy 

-5>6 

S6 

Q66- 

,£  , 

xy 

Q36- 

xy-' 


(29) 


where 


a' 

x 


l(k) 


x 


a*  }  =( 

y 

o'  ) 

xy 


y 

T 

V 


(k) 


fer. 


13 


-  ( Q23  r 


<5. 


36 


(k) 


r(k) 
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and 


Q..,QV 

Q..=Q..  =  Q..  --£  dJ  i  = 

i|  |i  i|  Q. 


‘33 


Q., 

Q..=Q  ..  =  7^ 

"  ''  Q33 


1,2,6  i  =  1 ,2,6 


1  =  1, 2,6  j  =  3 


» 


The  stress  resultants  and  stress  couples  are  obtained  by  integrating  the  stresses  over  the 
thickness  of  each  layer  and  summing  the  results  for  all  layers,  hence 

(k) 

'  dz 


N 

x 

N  h. 

a  „  k 

’  < 

a 

X 

•  N 

y 

-i  /  j 

a 

y 

N 

xyJ 

' 

k=i  hk_, 

a 

xyJ 

' 

M 

X 

N  h, 

_ _ _  *  k 

a 

X 

M 

y 

-I  / 

a 

y 

M 

xy 

k=i  hk_, 

a 

xy 

(30) 


(k) 


(31) 


Equations  (31)  therefore  represent  the  stress  couples  about  the  mid-surface  of  the  adhesive 
layer . 

3.4  Equilibrium  Equations 

Consider  a  thin  slice  cut  from  a  double  lap  joint  as  shown  in  detail  A  of  Figure  1 .  The 
surface  GH  for  this  case  is  not  a  free  surface  and  therefore  the  normal  stress  will  not  be 
zero  but  equal  to  a  .  The  shear  stress  along  GH,  however,  is  zero  because  of  symmetry. 
The  forces  and  moments  acting  on  the  upper  and  lower  laminate  of  such  a  slice  are  shown 
in  Figure  3.  Equilibrium  of  horizontal  forces  in  the  x-  and  y-directions  yields 


N  +  N  -  a°  =  0 
xu,x  xyu,y  xz 


N  .  +  N  .  +  o°  =  0 

xL,x  xyL,y  xz 


(32) 


and 
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(33) 


N  +N  -c  =0 
xyu.x  yu,y  yz 


N  .  +  N  .  +  a  =  0 

xyL-x  yL,y  yz 


Summing  vertical  forces  gives  the  relations  between  the  transverse  shear  forces  and  the 
normal  stresses  at  the  mid -surface  of  the  adhesive 


V  +  V  -o=0 

xu,x  yu,  y  z 


V  i  +  V  ,  +  a°  -  oc  =  0 

xL,x  yL,y  z  z 


Similarly  equilibrium  of  moments  in  the  x-  and  y-directions  gives 


M 

+  M 

-  V 

=  0 

xu,x 

xyu,y 

XU 

M  . 

+  M 

-v  , 

=  0 

XL,X 

xyL,y 

xL 

M  +  M  -V  = 0 
xyu,x  yu,y  yu 


M  ,  +M  .  -V  ,  =0 
xyL,x  yL,y  yL 


Equations  (32)  through  (36)  are  applicable  to  both  double  lap  and  single  lap  joints,  but  for 
the  latter,  a  in  the  second  of  equations  will  be  zero.  The  sum  of  the  stress  resultants  in 
the  upper  and  lower  laminate  represent  the  total  applied  loading  for  the  single  lap  joint 
and  one  half  the  applied  loading  for  the  double  lap  joint.  All  moments  are  again  calcu¬ 
lated  with  respect  to  the  adhesive  mid-surface.  It  must  also  be  remembered  that  in  a 
double  lap  joint  the  lower  laminate  represents  only  one  ha  if  of  the  laminated  plate  and 
that  the  transverse  shear  forces  and  moments  of  the  entire  plate  with  respect  to  its  middle 
surface  are  zero. 
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4.0  NARROW  UNIAXIALLY  LOADED  JOINTS 

4.1  General  Remarks 

Most  of  the  joints  tested  as  part  of  this  contract  were  1 .0  inch  wide  uniaxial  ly  loaded 

specimens  and  therefore  it  appeared  reasonable  to  assume  the  joints  to  be  in  a  state  of 

plane  sfress.  This  means  the  stresses  a  ,  a  ,  and  a  are  assumed  to  be  zero  while  the 
r  y  xy  yz 

remaining  stresses  are  constant  along  the  width  of  the  plate,,  i.e. 

.oW.0 

x,y  xz,y  z,y 

The  adhesive  material  will  be  assumed  isotropic  while  the  laminates  are  assumed  to  consist 
of  orthotropic  layers  (Section  3.1).  Half  of  the  adhesive  layer,  that  above  the  reference 
surface,  will  be  considered  as  part  of  the  upper  laminate  while  the  other  half  is  taken  as 
part  of  the  lower  laminate.  The  stacking  arrangement  of  the  laminated  plates  need  not  be 
symmetrica!  about  the  middle  surface.  The  joint  is  assumed  to  be  clamped  at  a  distance 
from  the  splice  plate.  This  assumption  does  not  affect  the  derivation  of  the  governing 
differential  equations,  of  course,  but  will  be  used  to  establish  the  necessary  boundary 
conditions. 


4.2  Single  Lap  Joints 

4.2.1  Effect  of  Transverse  Shear 

Because  of  the  plane  stress  assumptions  most  of  the  exp;essions  in  Section  3.2  are  simpli¬ 
fied  considerably.  For  a  laminate  consisting  of  a  single  layer,  equations  (15)  may  be 
written: 


axz  °  °xZ  ('  -  h^)+  t0®**  +  Q45V  f<2) 

V  =  °yz  ('  -  iT^)+  [Q45t(>x  +  Q44V  f(z> 

and  since  for  plane  stress  a  =  a°  =  0,  one  obtains: 

yz  yz 


(37) 
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<e^sassj3»jss mwavwwww 


O  /.  z  \  ,  1 


a  =  a  1  -  A-  hr^-cp  f(z) 

xZ  xZy  hN;  s55** 

Equation  (38)  may  be  generalized  for  a  multi-layer  laminate  by  writing  it  in  the  form: 


a(k)=o°  (l 
xz  xz  \ 


i^)+  ^kyf<z>  +  b2+°(k) 

L  55 


Continuity  of  shear  stress  at  the  interfaces  of  adjacent  layers  gives: 


,(k*»  -  <w  +  i  i  > 
1  +7kTTi?i)  f(hk> 

.  55  55  J 


and  as  before 


a(l)  =0 


Since  the  shear  stress  vanishes  at  z  =  h. one  finds: 

N 


b= 


The  transverse  shear  force,  V^,  is  obtained  by  integrating  the  shear  stress  over  the  thick¬ 
ness  of  the  laminate,  or: 

*VkzlV.+A55‘<’x  <42 

where: 


A55  ’  5bhN  +  2  ife  [F(hk>  -  F(hk-.»  +  °<k)<hk  *  hk-,> 

k=l  b55 

and  F(z)  is  as  defined  in  Section  3.2.  Solving  forcp^  from  equation  (42)  and  substituting 
into  (39)  yields: 


a(k)  =  G(k)c°  ±  G^V 
xz  1  xz  2  x 


in  which: 


'  -4  v 


tE  | 

£  6- 

I  r 

i  { 

gr  I 

I  I 

I  * 

&  l' 

I  «, 

&  * 

i  r 

i  *; 


5l«<k)=,-^-5hNG2(z)<k) 


r  M(k)_  1 

G2 z  "A 


55 


1 

55 


f(z)  +  bz  +  a^ 


(k)  (k) 

From  the  second  of  equations  (22),  setting  a'  "  0  and  substituting  for  a^z  from  equation 
(43),  one  has  for  the  transverse  shear  strain: 


e<k)  =  S(k)[G^k)a°  ±  G(k)V  ] 

xz  55  1  xz  2  x 


(44) 


f 

is» 


m 


m 


g<  f 


m  i 

I 


■i 


E 


Substituting  this  into  equation  (10)  and  integrating  with  respect  to  z  gives  the  displacement, 
u 

.  +Gik)°xz±G2>Vx  <45> 

o 

where 


u(k)  =  u°-  f  w,  dn+GV'V'  ±  GT'V 

J  x  1  xz  2  x 


G  (z)(k)  =  S<k>  f  G  (n)(k)dr,  +  c!k) 
1  ^  J  I  I 


.0) 


i  =  i,2 


(46) 


It  follows  from  (45)  that  the  constant  c;  y  must  be  zero.  Compatibility  of  displacement  at 
the  interfaces  requires: 


G  (hk)(k)  =  G  (hk)(k+1)  j  =  1,2 


(47) 


and  hence  the  remaining  constants  can  be  determined.  Differentiating  the  displacement 
(45)  with  respect  to  x  gives  the  longitudinal  strain: 

z 

-  f  w,  dr)  +  G(k)a°  ±  G^V  (48) 

J  xx  1  xz,x  2  x,x 


Jk)  -  ° 

e  -  e 
x  x 


— (k)  (k) 

It  must  be  noted  that  the  shear  functions,  G.  ’  and  G.  ,  do  not  have  the  same  meaning 

1  1 

os  in  Section  3.2. 
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4.2.2  Variation  of  Norma!  Displacement 

The  integral  in  equation  (48)  will  not  be  determined  consistent  with  the  assumed  transverse 
shear  deformation.  The  stress  equilibrium  equations  for  narrow  uniaxiaily  loaded  plates 


reduce  to: 


„»♦<,«  =0 
x,x  xz,z 

0<k>+c<k>  =0 

z,z  xz,x 


(49) 


(k) 

Integrating  the  last  one  with  respect  to  z  and  substituting  for  a  ,  gives  the  normal  stress 

XZ 


in  the  k^1  laver: 


7<k>=a°  +  P<kV>  ±P<k>V 

z  z  1  xz,x  2  x,: 


(50) 


where  a°  is  the  adhesive  normal  stress  at  the  reference  surface  and  the  normal  stress 
z 

functions  are  given  by: 


Pl(z)(  )  ~  2hN  *  Z  "  2  hN  *Vz) 


(k) 


P2(z) 


(k)  _  1 


55 


1  ri  \  .  1  I  2  (k)  (k) 

^)Hz)  2hz  +a  Z+P2 

55 


Continuity  of  normal  stress  requires: 


ww- w5*" 


(51) 


It  follows  from  equation  (50)  that  p„  =0.  The  remaining  pi  can  be  obtained  from  the 

(k)  1  (k)  1 

above.  Since  the  stresses  a  and  a  are  zerc.  one  may  write  for  the  normal  strain  in 

th  y  *y 

the  k  layer: 


„(k)  _  JO  _  s(k)  (k)  +  Ak)  flc) 

"z  "W'z"S13ax  +S33az 


(52) 


where: 


[S]  =  [Q]"1 
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Differentiating  with  respect  to  x  and  substituting  the  first  of  equations  (49)  gives: 


««-*«.«  +S«o« 

xz  13  xz,z  33  z,x 


(53) 


Substituting  equations  (43)  and  (50)  one  obtains  after  integrating: 


w 


lk)  =  w?  +pHf°  +Pj)a°  ±  P(k)V  iPjV  ±P(k)a° 

x  3  xz  4  xz,xx  5x6  x,xx  7  z,x 


(54) 


where: 


P3(z)(k>  =  S<.3^-5hNPyk>  +  P3k) 


P4(z)(k)  =  s<k)  j  P^nlcm  +  pJ* 


:(k) 


P5W(k>  =  -^ 

D  55 


1 


7M 

55 


f(z)  +  bz 


+Dw 


P6(d<k)  =  S33/P2W*'  +  p 

o 

P^W-sg^pW 

At  the  reference  surface  (z  =  0),  w^  =w°  ;  and  hence: 


00 

6 


,0) 


p:  -0  j  =  3,4, . . .7 


(k) 


The  remaining  constants  p.  are  determined  from  the  condition  that  the  slope  must  be  equal 
for  two  adjacent  laminae  at  their  interface,  or 


p  <hki(k)  =  p  (hk)(k+,) 


(55) 


Differentiating  equation  (54)  with  respect  to  x  and  integrating  with  respect  to  z  gives: 
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f  w  dri  =  zw°  +  p£k)a°  +  Pok)a°  ±  P?J?V  ±  P(k)V  ±  P(k)a° 
J  W/xx  'xx  8  xz,x  9  xz,xxx  10  x,x  1<  x,xxx  12  z. 


XX 


(56) 


where 


z  k-1  h 

P.+5(z)^  =  J  Pj(ri)^dr|+^  J  P j(r|/m^drj  j  =  3,4,  ...7 


'k-l 


m=l  h 


m- 


4.2.3  Laminate  Stress  Resultants  and  Couples 

Since  and  are  assumed  to  be  zero  and  only  the  longitudinal  strain,  e^.  \  is  of 
y  xy  * 

interest  here,  one  may  solve  for  the  in -plane  strains  from  equations  (29).  This  gives 


e 

x 

(k) 

r 

a 

x 

(k) 

f  ' 

z, 

(k) 

4 

•  e 

y 

.  =.[3|<k>. 

a 

y 

•  _  i 

Z2 

.  aw+. 

z 

T2 

€ 

xy- 

a 

xy 

•Z3' 

T3- 

(57) 


where: 


Q  = 


qh 

Q12 

Q16 

Q,„ 

Q_„ 

Qn( 

12 

22 

26 

.Q16 

^26 

lo 

V 

Z2 

>  =  |Q|' 

Z3J 

Q 


13 


fS. 


23 


‘36 


V 

V 

X 

T„ 

- 14 

a* 

2 

y 

T_ 

a1 

3 

xy' 

Frcm  the  first  of  equations  (57)  one  finds: 

<k)_  1  rJk)+7(k)  (k)  T(k), 

°x  "|(k)Ux  Z1  °z  "T1  ] 


11 


30 


and  after  substituting  (48),  (50),  and  (56),  this  becomes: 


a*k) "  sW  Ce°  -  -  T<ik>3 + .  +  rJV 


Q 


11 


1  J  1 


*z,x  "2  UXZ, 


XXX 


+  °  r  »(k)  o 

+  a  i-  R  'a  +  K  - 
3  z  A  z,xx*  5 


iR^’v  ±R<wv 

X/X  6  x,xxx 


(58) 


where 


R,(4W  -  CG,(2)(k)  -  PB(z/kl  +  z|k)Pt(2)(k)]  /gW 

R2(m)<k)  =  -P9<z)<k=)/g(k) 

R3<z>(k,  =  z<k>/5<k> 

R4(z)(ll)=-PI2(2)(k)/Q<k) 

R5(z)(k)  =  [G2(2)(k)  - P|0(2)(k>  +  2<k,P2(2)<k)]  /g(W 

Vz>W  =  -pntt<k,/5<k) 

For  the  case  under  consideration  the  equilibrium  equations  (32)  through  (36)  reduce  to: 


N  -  o°  *  0 

xu,x  XZ 

NxL  x  +  *°  =0 
XL,x  xz 


V  -  a°  =  0 . 

XU,X  2 


(59) 


xl,x 


+  a°  =  0 


(60) 


and 


M  -  V  =0 

xu,x  xu 

M  .  -  V  =  3 

xL,x  xL  J 


(61) 


Adding  the  first  and  second  of  equations  (59)  and..  (60;,  one  } 


has: 
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N  +  N  .  =  0 

xu,x  xL,x 


V  +  V  ,  =  0 

xu,x  xL,x 


and  therefore: 


N  +  N  .  =  Constant 

xu  XL 


V  +  V  .  =  Constant 
xu  xL 


The  first  constant  is,  of  course,  equal  to  the  applied  axial  load,  N^;  while  the  second 


constant  is  zero  from  considerations  of  symmetry,  so  that: 


V  =  -V 

XL  XU 


Also  by  adding  equations  (61)  and  substituting  the  above: 


M  +  M  “0 
xu,x  xL,x 


which  means' 


M  +  M  .  =  M  -  Constant 
xu  xL  x 


With  the  use  of  equations  (60)  and  (64),  the  longitudinal  stress  (58)  may  be  written  in  the 
general  form: 


M  _  1  r  °  o  T(k),  R(k)  o  R(k)  o 

Tx  -=TniC£x'ZW'xx"Tl  J  +  Rl  Jxz,x +  R2  axz,, 


+  [R^  +  R^V  +  [R*^  +  R^]V 

3  0  xu,x  4  6  xu,xxx 


which  is  applicable  ro  both  upper  and  lower  laminate.  When  the  above  is  substituted  into 
the  first  of  equations  (30)  and  (31),  one  obtains  the  stress  resultants  and  stress  couples  on 
the  lanrnatss: 


N_.uo  Du  o  nu  o  _  u  o 

-A.,e  -B..w,  +  R,.a  +  R.na 

XL*  11  X  11  XX  11  XZ,X  1 2  xz,xxx 


+  Z“v  +  Z1'  V  -  N' 

II  xu,x  12  xu, xxx  xu 


■ 


32 


where 


N  L  =  Anew  “  Bnw°  +Rna°  +  rLt° 
x*-  x  •'  xx  11  xz,x  12  xz,xxx 

+  zi’.v  +  Z^V  -N* 

11  xu,x  12  XU/XXX  xL 

+  Z21Vxu,x  +  Z22Vxu,xxx  'Mxu 

MxL  =  Bnex  '  Dn*°„  +  *i,o°  +  Rx,o° 

XL  1 1  x  1 1  xx  21  xz,x  22  xz,xxx 

+  Z21Vxu,x+Z22Vxu/xxx"Mk 


N 


An=±X  *W(hk- W 

k=l  Q11 

Bn=4l  =k)(hk~hk-i> 

k=l  U11 


N 


k=l  «»  1 


N 


r,, =±  y 

t  / 

k=i  h 


v(k) 


n  1  z,  f  W dz 


N 


*12  = 


‘2  / 

k=l  hk./ 


R2(x) 


(k)  dz 


R21 


N  \ 

1  [  «I«(fc).-. 

w  hk-/ 


(68 

cont'd) 
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N  nk 

**22  =  ±  ^  z 

k=1  h^ 

N  T(k) 

Ni=±  S  prV’k-i* 

k=i  Qn 


n  T(k) 

mx_±2Z  S[ky(hk'hk-1) 
k=l  Q11 


N  \ 


hr*l 


k=1  h, 


J  [R3(2){k;  +  R5(?.)(k)]c 


N  "k 

Z12  =  ±I  /  f- 

k=ih.  ; 


[R4(z)(k)+R6(z)(k)]dz 


-k 

[r<Jz)(k)  +  RJz)(k)]zdz 


"21  x  Z.  J  L  3  *'5V 

k=i  h,  ; 


N  k 

Z22  =  ±I  /  tV'2>(k)+R6(2)(l<>:lz')z 

k=i  h,  : 


In  the  above,  the  superscripts  u  and  L  are  emitted  for  convenience. 


4.2.4  Governing  Differential  Equations 

Differentiating  equations  (68)  and  substituting  into  the  first  of  equations  (62)  and  into 
equation  (65),  gives  in  matrix  form. 
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Bu  Dn 


11 

R12 

zn 

Z12 

21 

R22 

Z21 

Z22. 

xz,xx 

o 


where: 


A11  A11+A11 

Bn  =  Bn  +  B11 


Solving  for  e  and  w,  from  equation  (69)  yields: 

X*X  XXX 


XU,XXXXJ 


o 

e 

x,x 

R'n 

R12 

Z11 

Z1  2l 

o 

-w, 

XXX 

-R21 

*22 

Z21 

CM 

IN 

XU, XX 


in  which: 


R11  R12  Z11  Z12 
R21  R22  Z21  Z22 


An 

Bn' 

-1 

Bn 

°n. 

xu,xxxxl 


R11  R12  Z11  Z12 
R21  R22  Z21  Z22 


Substituting  the  first  and  third  of  equations  (68)  into  the  first  of  equations  (59)  and  (61), 
respectively,  and  using  equations  (70),  one  obtains  the  following  two  coupled  differential 
equations: 


where: 


o°  +R*  a°  +  R*  a°  +  Z?,  V  +Z*  V  =0 

xz  11  xz,xx  12  xz,xxxx  11  XU, XX  12  xu,xxxx 

R*  a°  +  R*  a°  +  V  +  Z*  V  +  Z*  V  =0 

21  xz,xx  22  xz,xxxx  xu  21  xu,xx  22  xu,xxxx 
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II 

*•— 

CxZ 

-RU  - 
R11 

-AU 

A11 

Rn 

-6n 

*21 

R*2" 

_  RU  . 
K12 

-Al! 

R12 

'Bn 

*22 

R*2.“ 

.  RU  . 
K21 

'Bn 

Rll 

-°n 

*21 

R*22- 

■R22  ' 

R12 

-Dn 

*22 

zn  = 

'  Zll 

-An 

2.1 

-Bii 

1  Z21 

z*  =  -  -  A  Z  -  B  Z 

12  12  11  12  11  22 

7*  =  7U  du  ~y  _  nu  "7 

a21  21  °11  ^11  U11  "21 

7*  =  7U  _  au  "7  _  nu  "7 

22  22  11  12  11  22 
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i^gpj^t^g^ viyrr-r.w,  TrzFZSZ?* 7Zxsr~V7  -"' 


ir^^-^r^VsTj- 


. ,  o  o  o  o 

V  =  t-  a  +  t .  a  +  t7  a  +  tQ  a 

xu  5  xz  6  xz,xx  7  xz,xxxx  8  xz,xxxxxx 


(74) 


in  which: 


*5  '  S3*2  /  S2 


f6  “  ^21  +  Z22  f2  +  s3  V  ^  s2 


\?  -  (R*2  +  Z*2  f3  +  s3  f4)  /  s2 


s2=1+s3  rl 


s  =  Z*  +  Z*  t 
3  21  22  T1 


^8  *  "  Z22  VS2 


Differentiating  equation  (74)  twice  and  substituting  it  together  with  (73)  into  the  first  of 
equations  (71),  one  obtains  the  following  eighth  order  differential  equation  in  0°^: 


C.  a°  +C0a°  -tC~a°  +  C .  a°  +  cr°  =0  (75) 

I  xz,xxxxxxxx  2  xz,xxxxxx  3  xz,xxxx  4  xz,xx  xz 


where 


|  f 

#■  * 

* 


C1  s4  f8 


C2  "  Z?2  *4  +  S4  V 


C3  ~  R12  +  Z12  f3  +  s4  *6 


C4  R11  +  Z12t2+S4t5 


s  =  Z*  +  Z*  t 
S4  11  12 'l 


The  general  solution  to  the  differential  equation  (75)  may  be  written  in  the  form: 

8 


xz 


-S>, 


X.x 


i=l 


(76) 
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(77) 


where  X.  are  the  roots  of  the  equation: 

c]  \8  +  c2\6  +  c3\4  +  c4\2+  1  =0 


By  letting  r\  -  X  and  dividing  through  by  C  ^  this  becomes 

ti4  +  ^  r?  +  C2  r}  +  C3  r\  +  =  0 

where 


i  =  1/2,3 


c  =  — 

4  C 


1 


(78) 


In  general  the  solution  of  equation  (78)  will  yield  a  pair  of  complex  conjugate,  or 
imaginary  roots,  and  hence  the  form  in  equation  (76)  does  not  lend  itself  for  numerical 
computations.  Considering,  for  example,  the  ccse  in  which  equation  (78)  has  two  real 
roots,  one  of  which  is  positive  and  the  other  negative,  and  a  pair  of  complex  conjugate 
roots,  the  solution,  equation  (76),  may  be  written  in  the  alternate  form: 

Q  a.jX  -a,jX 

a  =A,  e  +  A_  e  +  A_  cosB^x  +  A,  sinB^x 
xz  I  2  3  3  4  3 

V 

+  e  (A^  cosB^x  +  A^  sinB,-x) 

-oux 

+  e  ~  (Ay  cosB^x  +  Ag  sinB^x)  (79) 


where  a.  is  the  real  part  and  B-  the  imaginary  part  of  the  complex  root: 


X.  -  a.  +  '8- 

!  I  i 


In  general  one  may  write  equation  (76)  in  the  form 


o 

a 

xz 


(80) 
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'."here  each  function  cp.  has  one  of  the  following  five  forms: 


a.x 
Q  I 


COS0  .X 

sin8.x 

I 

a.x 

e  *  cosp.x 
a.x 

e  *  sinp.x 


For  the  case  when  the  joint  configuration  is  symmetrical  with  respect  to  the  middle  surface 
of  the  adhesive  layer,  i.e. 


Aii  =  An 

8u  =-Bi, 

Dn=Dii 


it  can  be  shown  that  Z^,  Z^,  R^]'  anc*  ^22  are  zero  anc*  *wo  differential  equations 
(71)  become  uncoupled.  Instead  of  equations  (80),  one  now  has: 


a°  =  Y  A.  cp.(x) 
xz  ZJ 

i=1 

8 

V  A.  cp.(x) 

XU  ZL  I  M 


(81) 


i=5 


In  addition,  one  finds  for  the  symmetrical  case: 


R11  R12  =  Z21  Z22  ° 


39 


4.2.5  Boundary  Conditions  and  Solution 

The  eight  constants,  A.  in  equations  (80)  or  (81),  must  be  determined  from  the  boundary 

I 

conditions  of  the  problem .  The  first  four  boundary  conditions  are  as  follows: 


OO  O  _  A 

a  =  a  (0)  =  0 

xz  xz 

a°^  =  a°  (L)  =  0 
xz  xz 


V°  =  V  (0)  =  0 

XU  XU 


(82) 


V  -V  (L)  =  0 

XU  XU 

In  other  words,  the  adhesive  shear  stress  and  the  transverse  shear  forces  vanish  at  x  =  0  and 
x  =  L,  In  the  Goland-Reissner  analysis,  the  adhesive  shear  stress  would  have  its  maximum 
va  ue  at  these  points,  which  is,  of  course,  a  violation  of  the  boundary  conditions. 

The  loads  and  moments  at  the  ends  will  provide  the  remaining  boundary  conditions.  It  will 
be  necessary,  however,  to  determine  the  strain  and  curvature  at  the  reference  surface 
first.  From  equations  (70),  after  integrating: 


e  =  R . .  a°  +  R10a°  +  Z..V  +  Z19V  +  Ac 

x  11  xz,x  12  xz,xxx  11  xu,x  12  xu,xxx  9 

-w,  =  R0,  a°  +  R„9  a°  +  Z_.  V  +  Z  _  V  +  Ain 

xx  21  xz,x  22  xz,xxx  21  xu,x  22  xu,xxx  10 

Substituting  this  into  equations  (68),  one  obtains: 


(83) 


N  = 
xu 


R*  a  -R*  a 


-  Z*  V 


-  Z*  V 


11  xz,x  12  XZ,XXX  11  XU,X  12  XU, XXX 


+  AU  A  +  BU  A  -  N ' 
A11A9  B11A10  N 


XU 


N=R?*a°  +  R?ia°  +  Z**  V  +Zt*V 

XL  11  XZ,X  12  XZ,XXX  11  XU,X  12  XU, XXX 
+  A11  A9+BUAIO'Nk 


(84) 
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With  the  additional  two  integration  constants,  and  Ajq,  there  are  ten  undetermined 
coefficients  and  therefore  six  more  conditions  are  required  for  a  solution.  Although  it 
appears  that  eight  more  boundary  conditions  (two  loads  and  two  moments  at  each  end)  are 
available,  only  six  of  them  are  independent.  Because  the  first  of  equations  (62)  and 
equation  (65)  were  used  in  the  derivation  of  the  governing  differential  equations  (71), 
only  three  load  and  three  moment  conditions  can  be  specified.  As  shown  previously: 
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M  +M  ,  =M 

XU  XL  X 


N  +  N  ,  =  N 
xu  xL  x 

where  N  is  equal  to  the  applied  load  and  therefore  a  known  quantity.  M  represents  the 

X  X 

total  moment  about  the  reference  surface  and  can  be  determined  from  the  condition  that 
the  slope  is  zero  at  the  ends  and  at  the  center  of  the  joint. 


In  the  region  AC  of  the  lower  laminate  (Figure  1),  the  strain  and  curvature  will  be 
assumed  constant  and  equal  to  and  w °  ,  respectively.  The  slope  at  C  now  becomes: 


-o  _  .  -o 
w,  =  L.  w, 
x  1  xx 


(85) 


The  load  and  moment  will  be  constant  also,  and  are  given  by: 


NxL'Nx'An"<-BU»L-NxL 

MxL-Mx  =  Birex-Dn»'xx-Mk 


(86) 


By  eliminating  e°  from  equations  (86)  a.  .  bstituting  (85),  one  obtains  the  relation: 


B 

M  +  M'  --il(N  +  N’.)  =  -dS\  w° 


11  -o 


xL  L 

A11 


xL 


11  "'xx  L 


■  w. 


1 


(87) 


where: 


°11  “Dll  A 


1J_ 

11 


In  a  similar  manner  one  finds  for  region  DE: 


-L  _  -l. 

w,  -  -  L_  w, 
x  2  'xx 


(88) 


and: 
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•.  i 


;  f 


r  t 


DU  f\U 

“11  L>1 1  — L 

M  +M'  — —  (N  +  N'  )  =-7— w, 

x  xu  .u  x  xu  L_  x 

A,,  2 


(89) 


An  expression  for  the  slope  in  region  CD  is  obtained  by  integrating  the  second  of  equations 
(83)  with  respect  to  x: 

w?  =  -  R01  a°  -  R_0  a°  -  Z„,  V  -  Z.0  V  -  A.-x  +  A*. .  (?C) 

x  21  xz,x  22  xz,xxx  l\  xu,x  22  xu,xxx  10  11 

Before  matching  slopes  on  both  sides  of  points  C  and  D,  an  inconsistency  of  minor 
importance  must  be  pointed  out.  According  to  equation  (54),  the  slope  varies  through  the 
thickness  of  the  laminate  in  region  CD.  However,  beam  theory  was  used  in  region  AC  of 
the  lower  laminate  and  in  region  DE  of  the  upper  laminate  to  derive  equations  (87)  and 
(89) .  It  would  therefore  be  better  to  use  the  average  slopes  of  the  laminates  at  the  ends  of 
region  CD  .  These  average  slopes  follow  from  equations  (54),  (60),  (64),  (82),  and  (90), 
hence: 


w 


x  <P4AV  ‘  R22)axz,xx  +  (P6AV  '  Z22)Vxu,xx  +  A11 
'x  =  (P4AV  '  R22)axz,xx  +  (?6A  V  "  Z22)Vxu,xx  "  A10  L  +  A'll 


-=?  ,o 


-L 

w 


(91) 


v.here: 


P4AV=ilP4<h,><,)  +  P4<hN><N,| 

P6AV  ‘I  I  W0>  1  +  P7<hl><'>  +  P7(hN><N>l 

Substituting  the  above  into  equations  (87)  and  (89),  one  obtains,  after  rearranging: 

<P4AV  "  y°xz,xx  +  (P6AV  '  Z2^Wxu,xx  +  A1 1  +  =TT  Mx  = 


11 


’I 


6lll 


m;l-t-<nx+n;l> 


n 


(92) 


and: 
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Equation  (82),  (92),  (93),  end  (94)  represent  a  system  of  twelve  conditions  with  twelve 
unknowns.  In  matrix  form  one  may  write: 


44 


i«> 


4<L> 


1,1  Hl,2 

-  -  -  -  -  Hl,12 

2,1  H2,2 

H2, 12 

V  (0) 

XU 


V  (L) 

XU 


N  (0) 

XU 


N  (L) 

XU 


N  +  N  . 

XU  XL 


M  (0) 

XU 


M  (L) 

XU 


eq.  (93) 
eq.  (92) 


M  ■+  M  . 

I  xu  xLj 


_H  12,1  H 1 2 , 2 


H12,12J  [AilJ 


The  conditions  involving  the  transverse  shear  force  are  obtained  from  (74)  qnd  those 
involving  the  stress  resultants  and  couples  are  given  by  equations  (84) .  Equations  (95)  may 
be  written  in  the  comoact  form:  , 


[H]  (A*}  =  (C*} 

where  the  coefficients  of  [H]  are  given  by: 


H,  .  =  cp.(0) 


-  cc.  (L) 
=  t:(0) 


*’iW- 


=  Jj(0) 
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H,  .=H.O) 

6/l  I 

H_  .  =  0 
7,1 

H8  ;  =  P;<°> 

0,1  I 

H  =-  p  (L) 

7,  l  I 

H10,i  =  ^P4AV  " *22^fpj +  ^P6AV  "  Z22^f^ 

HH,i  =  <AV-R_22>’1i'(0>  +  (PMV-Z22>*i'(0) 

Vi'° 

*i  “  Vi +  Vi' +  ,7tpiv + vj1 

pi =  -  R2icpi "  «*a*r  “  Z2i  ■  i "  “22  ^i” 


H.  . 

M 

=  0  i  =  1,2, 3, 4 

1  =  9,10,11,12 

Hi,9 

n 

X 

> 

M 

II 

O 

i  =  5,6,7,8 

H-5, 10 

'H6,io-An 

H5/!i 

=  H6.11  =  Bli 

\lQ 

=  A1I 

H7,ll 

H8. 10 

"H9,I0  =  S11 

r*~ 

00 

X 

■"9.11 

H9,9 

=  -1 

HS,12 

H10,9 

H10,10 
H10,11 
H 1 0, 12 


H1 1 ,9 

Hn .  io 

Vll 
H1 1 , 1 2 
H12,9 
H 1 2, 10 
H12, 1 1 
H12, 12 


=  0 

=  0 

=  1 
» -1 


-  0 


and  the  coefficients  of  the  column  matrix  {C*}: 


cr 

I 


=  0  i  =  1 ,2,3,4 


=  N' 
xu 


-  N  +  N' 
x  xu 


=  N  +  N'  +  N' . 
x  xu  XL 
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C8  =C9  =  ML 


4  B11 

C*  =—  M‘  — — ■  (N  +N'  ) 
10  ~u  xu  ,u  x  xu 
D,,  i  A., 


L1  I  Bu 

’■=i~  Mk--r<N*  +  Nk> 
Dn  L  An 


C*  =M*  +M' 

12  XU  XL 


Solving  the  set  of  simultaneous  equations  (96)  yields  the  unknown  coefficients  A?  and 

therefore  A.  ano  M  . 

i  x 

{A*}  -  [H]“ 1  [C*}  (97) 

After  determination  of  the  coefflcieni's  A.,  ihe  sheer  and  normal  stresses  in  the  adhesive 

I 

may  be  calculated  from  equations  (80)  and  (60),  or  in  the  symmetrical  case  from  equations 
(81)  and  (60).  The  interlaminar  shear  and  normal  stresses  in  the  adherends  nay  be 
obtained  from  equations  (43)  and  (i>G),  and  the  longitudinal  stresses  in  the  lanr.incte  are 
given  by  equation  (67). 

4.3  Double  Lop  Joints 


For  the  upper  laminate,  the  cna lysis  is  identical  to  that  for  the  single  lap  joint,  but  some 
minor  modifications  ere  required  with  respect  to  the  lower  laminate.  As  a  result  of  sym¬ 
metry,  the  laminated  plate  os  a  whole  w;i!  not  incur  any  bending  and  therefore  it  will  be 

reasonable  to  assume  the  sheaf  function,  co  ,  to  be  zero  for  the  lower  leminate.  With  this 

x 

assumption,  the  transverse  sheer  stress  will  vary  linearly  through  ihe  laminate  according  to 


the  relation: 


0(k)=a°  [1-4 

xz  xz  \  h. 


(k) 

while  the  shear  strain  s'  '  becomes: 

xz 
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,«-S«  1  -,i-)o0 

xz  55  y  h^y  xz 

Substituting  this  into  equation  (10)  and  integrating  with  respect  to  z  yields  for  the 
longitudinal  displacement: 


(k)  _  o  r  _(k)  c 

u  -u  -  /  w,  dn  +  G.a 
J  x  '  1  > 


where 

G,w  -s55 

The  longitudinal  strain  is  obtained  by  differentiating  the  displacement  (100)  with  respect 
tox.  Hence: 


(k)  _  o  r  ,  .  _(k)  o 

e  -  e  -  /  w,  drji-G.a 

X  X  J  XX  1  XZ,) 


The  shear  force  acting  on  the  lower  laminate  becomes: 

VxL="20xzhN  (102) 

The  normal  stress  in  the  k*^1  layer  is  found  by  integrating  the  second  of  equations  (49)  with 
respect  to  z  and  substituting  (98): 

=0°+<4,*(-2+5b)  (,03) 


Equation  (53)  is  valid.  By  substituting  the  expressions  (98)  and  (103)  and  integrating  the 
result  with  respect  to  z,  one  has  for  the  slope: 


,(k)  =  w?  +P«00  +  P<k>a°  ♦  P<k>a° 
x  x  3  xz  4  xz,xx  7  z,) 


where 


N 
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P4<2><k)  =  S33(-T+6^)+Pr 

p7(z)(k)  =  S<k)  z  ♦  P(7k) 

The  constants  pt\  p^,  and  p^  are  determined  from  conditions  of  slope  compatibility  of 

O  / 

adjacent  layers  at  their  interfaces.  From  the  above  one  obtains: 

Z 

/.  o  .  D(k)  o  D(k)  o  .  D(k)  o  /iak\ 

w,  dri  =  zw,  +  PQ  a  +  PI  a  T  Plo0  (1 05) 

'xx  xx  8  xz,x  9  xz,xxx  12  z,xx 


in  which,  as  before: 


p.+5(z)(k)  =  j  P.w00*,  +  2 


m=l  h 


n 

/m 

Pi 


wH* 


i  =  3,4,7 

The  equilibrium  equations  for  the  double  lap  joint  are  the  same  as  those  for  the  single  lap 
joint  except  that  in  place  of  the  second  of  equations  (60),  one  has: 


V  +  a°  -  aC  =  0 
xL,x  z  z 

The  longitudinal  stress  for  the  lower  laminate  may  now  be  written: 


(k)  _  1  I  (k)  _ 


?s(k)  lx 


£  -  ZW 


?  -T(k)  +  R^k>0°  +  R<k)a° 

xx  1  l  xz,x  2  xz,) 


+  r!)<)V  +R*k)V  nn7, 

3  xu,x  4  xu,  xxx  (10') 

(k)  (k)  (k) 

which  is  identical  to  equation  (67),  if  the  quantities  Pj  ,  R^  /  and  in  that  equation 
are  set  equal  to  zero.  Equations  (64)  and  (65)  are  no  longer  valid  but  the  first  of  equations 
(62)  still  applies  end  hence: 

N  +  N  .  =  N  (108) 

XU  XL  X 
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As  mentioned  previously,  however,  in  this  case  represents  only  half  of  the  applied 
load.  The  expressions  for  the  stress  resultants  and  stress  couples  given  by  equations  (68) 
remain  valid,  but  the  coefficients,  R..  and  zh,  must  be  modified  in  accordance  with 

l|  !| 

(107).  Since  the  slope  of  the  lower  laminate  is  zero  for  z  =  h^,  one  has  from  (104),  and 
the  first  of  equations  (60): 

w?  =  -  ?lJhJN)a°  -  P^(hw)(N)cj0  -P^(S  )(N)V  (109 

x  3  N  xz  4  N  xz,xx  7  N  xu,xx 


Differentiating  the  above,  substituting  into  the  first  and  second  of  equations  (68),  and 
using  (108)  yields: 


A  e°  + 
1 1  x 


R,1  +BllP3(hN)lN)|°xZ,x+  lR12+  B,1^V(N)1°°z,> 


+  Z..V  +  Z...+ B..  P_(hM)'”'  V  =  N  +  N’  +N*. 
11  xu,x  l  12  11  7  N  xu,xxx  x  xu  xL 


and  hence: 


e°  =  R11a°  +R10a°  +  Z..V  +  Z.„V 

x  11  xz,x  12  xz,xxx  11  xu,x  12  xu,xxx 


where 


+  a~  (N  +  N*  +  N' .) 
A  j  j  x  xu  xL 


Rn  = 


‘  Ki  +  BnP3(hN)*  *1  /Au 

■|R12  +  BllVhN^  *1  /AU 


zn  “ 


-Zi/An 


Z12  = 


Z12  + 


b„p^n)(n)| /A„ 


After  substituting  (109)  and  (1 10)  into  the  first  and  third  of  equations  (68),  and  using  the 

equilibrium  equations  (59)  and  (61),  one  obtains  two  fourth  order  differential  equations 

identical  to  (71).  The  coefficients  R*  and  Z*.  for  this  case  are  given  by: 

MM 
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If'  _ _ _ _ _ 


r  ? 


Rl,  “-Rn-AnR,i-Bnp3<hN><N) 

R12  =  'R12“AllRt2'Bll  PWN> 


U  _  DU  D  _  r\U  P^-/L  \(N) 

21  BnRn  Dn  P3(V 


R21  =  " R 


R22  '  R22  ‘  B11  R12  "Dll  P4®v/  ' 


(N) 


ZU  =  "Z11  ‘All  z„  -Bi,  P^N> 


Zi2=-Z12-AH  Z12 

Z2I  '  'Z21  “Bll  Z11  ‘Dll  P>N>(N) 


7*  =  _  7U  _  QU  7 

*“22  22  D11  H2 


The  governing  differential  equation  (75)  remains  unchanged  but  the  modified  definitions 

for  R*  and  Z*  must  be  used  to  determine  the  coefficients  C,  through  C  ..  The  solution  is 
1 1  1 1  I  4 

then  again  given  by  equation  (80). 

Equations  (90)  and  (109)  represent  the  slope  at  the  reference  surface  for  the  single  and 
double  lap  joint,  respectively.  The  two  expre3sions  become  identical  by  setting: 

R  =  PLfh  }<N> 

R21  P3(W 

R  -  PLfh  ){N) 

/2  P4  ”N 

7  =PLfh  ><N> 

Z22  /  N 

Z  =  A  =A'  =0 
21  M10  11  ^ 

for  the  double  lap  joint.  Similarly,  by  letting: 

1 


'll 


(N  +N*  +  N! . )  =  A 

”  “■  XL 


XU 


(111) 
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in  equation  (110),  the  expression  for  the  strain  e°  becomes  identical  to  that  given  by  (83) 
for  the  single  lap  joint.  The  correctness  of  (1 1 1)  can  easily  be  verified  by  substituting 
Ajq  =  0  in  the  seventh  of  equations  (95) .  Since  equation  (66)  is  not  applicable  and  the 
condition  of  zero  slope  at  the  ends  of  the  joint  (point  A)  is  satisfied  automatically,  the 
last  two  of  equations  (95)  are  discarded  for  the  double  lap  joint.  By  dropping  the  last 
two  columns  of  the  coefficient  matrix  [H]  (those  corresponding  to  A^q  and  A  j  j)  a  system 
of  ten  equations  with  ten  unknowns  remains;  the  solution  of  which  is  given  by  (97). 

4.4  Numerical  Results 

A  number  of  runs  were  made  w  rh  the  computer  program  BONJO  I  in  order  to  verify  the 
numerical  results  and  to  study  the  effect  of  certain  parameters.  No  experimental  data 
could  be  found  that  provided  information  with  regard  to  the  stress  distributions  within  a 
joint  as  virtually  all  such  data  deals  with  overall  joint  strength  or  average  shear  stresses. 
Comparisons  were  made,  however,  with  the  Goland-Reissner  analysis  and  with  finite 
element  analyses  performed  at  Lockheed.  The  results  are  presented  and  discussed  below. 

4.4.1  Adhesive  Shear  and  Normal  Stresses 

Adhesive  shear  and  normal  stresses  were  calculated  for  a  single  lap  joint  with  aluminum 
adherends,  using  the  Goland-Reissner  procedure  for  joints  with  relatively  inflexible 
cement  layers  as  reported  in  reference  1 .  The  same  joint  was  then  analyzed  in  detail  with 
BONJO  I  (symmetrical  version)  in  order  to  make  a  direct  comparison  between  the  two 
procedures  and  study  the  effects  of  transverse  shear  and  normal  stresses  in  the  adherends  on 
the  adhesive  stress  distributions.  These  effects  are  ignored  in  reference  1 .  Curves  (T) 
and  (2)  in  F  igures  4  and  5  show  a  direct  comparison  of  the  adhesive  shear  and  normal 
stresses  between  the  Goland-Reissner  procedure  and  BONJO  I. 

It  can  be  seen  from  Figure  4  that  the  peak  shear  stresses  do  not  occur  at  the  ends  of  the 
joint  and  that  the  magnitude  of  these  stresses  *s  considerably  lower  than  those  calculated 
with  the  Goland-Reissner  approach.  This  is  attributed  to  the  fact  that  the  transverse  shear 


53 


_ I _ 

TU’^T""  ?/ 1 


|(l)  Goland-Reissner  Analysis  (infinitely  rigid  adherends). 

-©  Lockheed  bonded  Joint  Analysis  (BONJO  I),  including  effect  — 
of  transverse  shear  and  normal  stress  ir.  adherends. 

Qn  BONJO  I,  neglecting  effects  of  transverse  shear  ard  normal 
stress. 

-©  30NJC  I,  neglecting  effect  of  normal  stress  but  accounting 
for  transverse  shear. 

( BONJOoI^ne|lecting  effect  of  transverse  shear  but  accounting  ' 


- *—  x(m.) 

FIGURE  **  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  IN  SINGLE  LAP  JOINT, 
COMPARISON  WITH  GOLAND-RE ISSUER  ANALYSIS 
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ADHESIVE  NORMAL  STRESS 


FIGURE  5  ADHESIVE  NORMAL  STRESS  DISTRIBUTION  IN  SEiGiB  LAI  JOINT, 
COMPARISON  WITH  GOLAND-RE ISSNER  ANALYSIS 


moduli  as  well  as  the  elastic  moduli  through  the  thickness  of  the  joint  have  finite  values 
in  the  Lockheed  analysis.  Curve  ^3)  shows  the  adhesive  shear  stresses  obtained  with 
BON  JO  I  if  both  the  transverse  shear  modulus  and  the  modulus  of  elasticity  of  the 
adherends  in  the  direction  normal  to  the  plane  of  the  joint  are  made  very  large  and  it  is 

shown  that  the  stresses  in  this  case  agree  much  closer  with  those  of  reference  1 .  Curve 

.  /■ 

represents  the  case  when  only  the  effect  of  transverse  shear  is  accounted  for,  while  curve 
(5)  takes  only  the  normal  stress  into  account.  A  somewhat  surprising  result  is  that  the 
larger  (xirt  of  the  difference  between  the  two  procedures  lies  in  the  fact  that,  in  the 
Goland-Reissner  analysis,  the  adherends  are  assumed  to  be  incompressible  in  the  direction 
normat  to  the  plane  of  the  joint.  A  comparison  between  the  two  methods  with  regard  to 
the  adhesive  normal  stresses  (peel  stresses)  is  shown  in  Figure  5.  Again,  most  of  the  dif¬ 
ference  is  caused  by  assuming  that  the  adherends  are  incompressible  normal  to  the  plane  of 
loading . 

The  effects  of  transverse  shear  and  normal  stresses  in  the  adherends  on  the  adhesive  shear 

stresses  have  also  been  investigated  for  an  aluminum  double  lap  joint.  The  joint  was  first 

analyzed  by  using  the  actual  isotropic  material  properties  of  aluminum.  The  resulting 

shear  stresses  are  shown  in  Curves  (?)  of  Figures  6  and  7.  These  curves  therefore 

represent  the  actual  shear  stresses  for  the  joint  under  consideration.  In  order  to  determine 

the  effect  of  transverse  shear,  the  same  joint  was  analyzed  again  after  arbitrarily  decreasing 

and  increasing  the  transverse  shear  modulus  G'  by  a  factor  of  10.  This  is  shown  in  Figure 

6  by  Curve;  (T)  and  (3)  respectively.  No  appreciable  change  in  shear  stresses 

resulted  when  the  modulus  was  increased  but  when  the  modulus  was  ducreased  from 
6  5 

4.0  x  10  to  4.0  x  10  psi  the  peak  shear  stresses  were  reduced  significantly.  The  effect 
of  the  normal  stresses  in  the  adherends  on  the  adhesive  shear  stresses  was  determined  in  the 
same  manner  and  the  results  are  shown  in  Figure  7. 

The  effect  of  using  finite  values  for  the  transverse  shear  and  normal  stiffnesses  in  the 
adherends  on  the  adhesive  shear  stress  distribution  in  a  composite  double  lap  joint  is  shown 
in  Figure  8.  In  this  case  BONJO  I  is  compared  with  the  modified  Goiand-fieissner  method 
described  in  Chapter  2.0.  The  magnitude  of  the  peak  shear  stresses  computed  with  the 
latter  method  was  more  than  50  percent  greater  than  those  calculated  with  BONJO  I. 
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FIGURE  7  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  IN  DOUBIZ  LAP  JOINT- 
EFFECT  OF  MODULUS  OF  EIASTICITY  NORMAL  TO  PLANE  OF  JOINT 
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FIGURE  8  EFFECT  OF  TRANSVERSE  SHEAR  AND  NORMAL  STRESS -COMPOS ITE 
DOUBLE  LAP  JOINT 


W  -  .  YVlt  WA  N  ^  _ 


These  results  substantiate  the  statement  made  earlier  that  an  analysis  based  on  the 
Kirchhoff  assumptions  could  lead  to  unacceptable  errors. 

The  presence  of  residual  thermal  stresses  may  become  important  in  the  analysis  of  composite 
joints.  These  stresses  are  caused  by  bonding  at  elevated  temperature  and  subsequent 
cooling  to  operating  temperature.  The  adhesive  shear  and  normal  stresses  in  a  boron/ 
aluminum  double  lap  joint  caused  by  a  temperature  differential  of  -175°f  is  shown  in 
Figure  9.  Since  the  thermal  shear  stresses  are  positive  on  one  side  of  the  joint  and  nega¬ 
tive  on  the  ether,  they  will  usually  cause  an  increase  of  the  maximum  adhesive  shear 
stress  in  a  loaded  joint.  This  is  illustrated  in  Figure  10,  which  shows  the  joint  subjected 
to  a  load  of  2250  lbs.  ,The  adhesive  shear  stresses,  due  this  load  ere  compared  with  those 
resulting  from  the  same  load  plus  a  temperature  differential  of  -175°F.  The  maximum 
adhesive  shear  stress  was  increased  from  3550  to  5400  psi  as  a  result  of  thermal  effects. 

Comparisons  with  finite  element  analyses  (FAMAS*)  were  made  for  a  few  lap  joint  con¬ 
figurations  and  good  agreement  with  BONJO  I  was  obtained  in  all  sases.  One  of  these 
cases  will  be  presented  here.  The  joint  configuration  and  the  material  properties  used  are 
given  in  Figure  1 1 .  The  joint  was  modeled  with  triangular  anisotropic  constant  stress 
elements  in  the  x-z  plane  as  shown  in  Figure  12.  Element  stiffness  properties  corre¬ 
sponding  to  a  plane  stress  condition  (a^  =  a ^  =  a  =  0)  were  used  in  the  program.  A 
comparison  of  the  adhesive  shear  and  normal  stresses  obtained  with  the  two  methods  is 
given  in  Figure  13.  Agreement  between  the  two  procedures  was  extremely  good  except  in 
the  immediate  vicinity  of  the  joint  edges  where,  as  expected,  it  was  not  possible  to 
approach  a  zero  shear  stress  condition  with  the  finite  element  analysis. 

4.4.2  Adherend  Stresses 

Axial  stresses  in  the  boron  laminate  and  the  titanium  splice  plate  for  the  joint  configura¬ 
tion  shown  in  Figure  11  were  calculated  with  BONJO  I  and  FAMAS  for  the  purpo»e  of 
comparing  results.  The  titanium  was  divided  into  four  equal  slices  and  the  stresses  were 


*Flutter  And  Matrix  Algebra  System;  Lockheed  developed  anclysis  programs. 
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FIGURE  9  ADHESIVE  STRESSES  CAUSED  BY  TEMPERATURE  DIFFERENTIAL  OF  -175°F 
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FIGURE  11  COMPOSITE  DOUBLE  LAP  JOINT 
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FIGURE  13  SHEAR  AND  NORMAL  STRESS  DISTRIBUTION  AT  CENTER 
LINE  OF  ADHESIVE  IN  COMPOSITE  DOUBIE  IAP  JOINT 
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calculated  at  the  five  locations  shown  in  Figure  14.  Finite  element  results  are  shown  at 
the  outer  surface,  the  midsurface  and  at  the  titanium/adhesive  interface.  The  stresses 
across  the  thickness  of  the  splice  plate  are  essentially  uniform  for  a  large  portion  of  the 
joint.  In  areas  of  appreciable  bending  BONJO  I  shows  generally  higher  gradients  than 
FAMAS,  but  it  appears  that  this  is  a  result  of  the  fact  that  the  FAMAS  model  was  not 
detailed  enough  in  these  areas.  Axial  stresses  in  the  boron  laminate  were  calculated  and 
plotted  in  Figure  15  for  the  three  locations  shown.  Good  agreement  with  the  finite  ele¬ 
ment  analysis  was  obtained.  Finally  Figure  16  shows  a  plot  of  the  transverse  shear  stress 
and  normal  stress  at  the  joint  cross  section  x  =  0.025"  which  is  near  the  point  of  maximum 
adhesive  shear  stress.  Finite  element  results  are  not  shown  here  but  agreement  with  FAMAS 
results  was  excellent. 


4.5  Joints  with  Non-Linear  Adhesive  Stress-Strain  Behavior 


4.5.1  General 


A  major  drawback  of  the  method  of  analysis  described  in  Sections  4.2  and  4.3  is  the  fact 
that  it  is  based  on  the  assumption  that  all  stresses  in  the  joint  remain  in  the  elastic  range. 
This  assumption  may  usually  be  justified  with  regard  to  the  stresses  in  the  adherend 
materials,  but  the  joint  configuration  is  generally  such  that  the  adhesive  stresses  will 
reach  the  proportional  limit  of  the  material  at  an  early  stage  of  loading.  For  this  reason  it 
became  necessary  to  extend  the  linear  analysis,  to  permit  non-linear  stress-strain  behavior 
in  the  adhesive  material.  Several  iterative  type  methods  were  investigated  to  solve  this 
complex  problem.  Some  of  these  are  discussed  below. 


4.5.2  Secant  Modulus  Approach 


The  portion  of  the  joint  defined  by  0  s  x  s  L  (see  Figure  1)  was  divided  into  regions  and 
secant  moduli  were  calculated  for  the  adhesive  material  based  on  the  average  stresses  in 
each  region.  These  secant  moduli  were  then  used  in  place  of  the  elastic  moduli  of  the 
adhesive  to  obtain  the  coefficients  of  the  governing  differential  equations  for  each  region. 
Conditions  of  continuity  and  compatibility  were  satisfied  between  adjacent  regions,  which 
together  with  the  boundary  conditions  at  the  ends  of  the  joint  were  necessary  to  obtain  a 
solution.  This  procedure  was  used  to  establish  basic  joint  design  parameters  early  in  the 


66 


Cloeed  Ft.rm  Analysis  at  Location  (^3^ 

Closed  r'o'ii  Analysis  at  Indicated  Locations 

Finite  Element  Strees  at  LocationA) 

Finite  Element  Stress  at  Locations) 

Finite  Element  Stress  at  Locaticr/V) 


AXIAL  STRESS  -  K3I 


Titanium  Splice  Plate 


<£•  Adhesive 


Li... 


Joint 


!  ' - -  Boron  Laminate 

ELASTIC  ANALYSIS 

Closed  form  Analysis  at  Locations  Indicated 
Finite  Element  Stress  at  Location  (T) 

Finite  Element  Stress  at  Location 
Finite  Element  Stress  at  Location  (3) 


Note:  Stresses  at  Location  (T) 
must  be  divided  by  10 


x  (in.) 


.10  .20 


•30  .40  .50 


.60  .70 


FIGURE  15  AXIAL  STRESSES  IN  BORON  LAMINATE, 
COMPOSITE  DOUBLE  IAP  JOINT 
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FIGURE  16  TRANSVERSE  SHEAR  AND  NORMAL  STRESS 
COMPOSITE  DOUBLE  LAP  JOINT 
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program  but  it  turned  out  to  be  unsatisfactory  for  joints  in  which  a  relatively  large  portion 
of  the  bond  layer  is  in  the  plastic  range.  Computational  difficulties  arise  when  matching 
boundary  conditions  between  adjacent  regions  because  of  the  discontinuous  nature  of  the 
analysis.  For  this  reason  development  work  using  this  approach  was  discontinued. 


4.5.3  Analysis  Based  on  Deformation  Theory  of  Plasticity 

The  theoretical  work  for  this  procedure  has  been  completed  and  is  described  below. 


Assuming  that  the  only  non -zero  stresses  in  the  adhesive  are  the  shear  stress,  a  ,  and  the 

xz 

normal  stress,  a^,  one  may  write  the  stress-strain  relations  in  the  adhesive  as: 
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where  e^,  e  ,  e^,  and  are  the  plastic  strain  components.  The  equivalent  stress. 


equivalent  total  strain  and  total  plastic  strain  are  defined  by  the  expressions: 
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which  are  related  by  the  equation: 
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With  the  assumption  that  the  volume  remains  constant  in  the  plastic  range  one  obtains 
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It  can  be  shown  that  the  individual  plastic  strain  components  are  given  by: 
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By  using  a  Ramberg -Osgood  representation  of  the  stress-strain  curve  in  the  yield  region 
one  has 
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and  since  n  and  Cq  _are  constants  for  a  given  material,  the  plastic  strain  components 

(116)  may  be  determined  if  the  stresses  a  and  a  are  known . 
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If  in  the  linear  analysis  procedure,  the  elastic  stress -strain  relations  a-e  replaced  by 
equations  (112)  and  if  the  plastic  strains  are  assumed  constant  through  the  .r.ickness  of  the 
adhesive,  the  two  governing  differential  equations  are  ro  longer  homogeneous  but  in  the 
form: 
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As  was  done  in  the  linear  procedure,  the  two  differential  equations  (118)  may  be  combined 
into  a  sit.gle  eighth  order  differential  equation: 

Z.  d!F. 

°  j.  r  °  ,  r-  °  ,  r-  °  ,  O  ST  i 

C,  a  +  C0  a  +Co0  +  C.a  +c  =  >  — - 

!  xz,xxxxxxxx  2  xz,xxxxxx  3  xz,xxxx  4  xz,xx  xz  Z_f  .  i 


where 


F0  =  Tj-e 
2  '1  xzp 

V’Vzp 

*4  =  ^21^1  "  ZllV5xzp 
F5  =  -  (Z^q2  -  Z^ri^p 


•_o  dx 
i-2 


(119) 
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P6  =  ^Z22T11  “Z12T13)exzp 

^7=  '  *Z22n2  "  Z12Vezp  . 

In  deriving  the  differential  equation  (119),  the  following  expression  was  obtained  for  the 
transverse  shear  force  in  place  of  (74) 

I 

V  =  tjj°  +  t  ,a°  +  t  -p°  +  iQa 

XU  DXZ  6  XZ,XX  /  XZ,XXXX  8  xz,xxxxxx 

q_e  +  e  +  r\-£.  + ‘Hoe  (120) 

5  xzp,xx  6  zp,xxx  /  xzp,xxxx  o  zp,xxxxx 

where 


VrV1^ 

4 


^6  s ,  ^2  "  n4 
4 

T17="^<,1lZ22"n3Zl2) 


4  • 


and  s^  and  s^  are  as  defined  in  Section  4.2.4. 


Since  a  direct  solution  of  equations  (120)  is  not  possible,  an  iterative  procedure  must  be 
used  so  that  the  piastic  strains  and  their  derivatives  will  be  treated  not  as  unknowns,  but 
as  known  quantities  obtained  from  the  previous  iteration.  Attempts  to  represent  the  plastic 
strains  by  continuous  functions  such  as  truncated  Fourier  series  or  power  series  were  unsuc¬ 
cessful  because  of  the  extremely  high  gradients  that  exist  near  the  edges  of  the  joint. 
Other  possibilities  were  investigated,  the  most  promising  of  which  appeared  to  be  the  use 
of  Green's  functions  coupled  with  a  numerical  integration  procedure.  Using  the  latter 
approach,  one  may  write  the  solution  of  the  differential  equation  (119)  in  the  form: 


a  x  / 

°°*z  ’  I  Ai  ”iw  +  /  I  F;<! 


o  5=2 


p.(s) d  G(V-5)  dg 

'  dx' 
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where  G(x  -  §)  is  the  Green's  function  ot  the  differential  operator: 


.8  A6  A4  J2 

^  '  =  C1  Tli +  C2  +  C3  7~4  +  C4  T"2  + 
dx  dx  dx  dx 


The  first  sumrpation  in  equation  (121)  represents  the  elastic  solution  to  the  problem.  The 
•form  of  the  Green's  function  is  dependent  on  the  type  of  roots  resulting  from  the  linear 
solution.  In  general,  one  has: 

4 

G(x  -§)=y  K  cp*(x-S)  (122) 

1  ,  n  n 

n=1 


in  which  each  of  the  functions  cp*  has  one  of  rhe  following  forms: 

1)  sinha(x 

2)  sina(x  -  §) 

3)  sinha(x  -  §)  cosp(x  -  §) 

4)  cosha(x  -  5)  sinP(x  -  §) 

corresponding  tc  a  real  positive,  a  real  negative,  and  a  pair  of  complex  conjugate  (or 
imaginary)  roots  of  equations  (78),  respectively.  For  convenience  it  will  be  assumed  that 
the  Green's  function  contains  one  of  each  of  these  forms,  hence: 

G(x  -  §)  =  Kj  sinha.j(x  -  %)  +  Kj  sinP^x  -  §)  +  sinha^x  -  §)  cosp^(x  -  §) 

+  cosha^(x  -  §)  sinp^fx  -  §)  (1 23) 


The  constants  K  are  determined  from  the  conditions: 
n 

dG(0)  _  d3G(0)  _  d5G(0)  _  „ 

r~  « 


dx 


dx' 


dx" 


d  G(0)  _  1 
dx7  'ci 

Higher  derivatives  of  the  Green's  function  are  given  by: 


(124) 


n  —  I 

^ a"  cosha.j(x  -  %)  +  (-1)  ^  cos^Cx  -  ?) 

dx 

+  (a  K_  +  b  K  .)  cosha,-(x  -  §)  cos0.(x  -  §) 
n  0  n  4  5  0 

+  (a  K  .  -  b  K-)  sinha^x  -  §)  sin0  _(x  -  §) 
n  4  n  «  j  0 


— ~  =  ^  1  ^  s'nbaj(x  "  §)  +  (-1)^  ^2*3  5'n^3^x  “  ?) 


m 


dx 


+  (a  K*  +  b  K  .)  coshctj-Cx  -  §)  sin0,-{x  -  ?) 
m  o  m  4  o  o 

+  (b  K0-a  K  .)  sinhac(x  -  §)  cos 0c(x  -  §) 
m  o  m  4  0  0 


025) 


where  m  is  an  even  and  n  is  an  odd  number,  and: 


a  =0 
.  o 


a  =  0_a  ,  +  cub  , 

n  5  n-i  5  n-l 


a  =  “  (3r°  1  +  a,-b  ,) 

m  0  m-l  o  m-l 


b  =  1 

o 


b  -  0  b  ,  -  cua  , 
n  5  n-i  5  r.-l 


b  "  l  “  S«;b  1 
m  0  m-l  o  m-l 


Equations  (125)  may  be  wr  itten  in  alternate  form: 


dnG(x  - 1)  _  V  „  ^  V  .  „ 

— — - S,nfe)e  =<”85*  +  g^CJe  sinSjX 

dx 


-CX5X 


-a5x 


+  g3n(?)e  COS05X  +  g  fe)e  5^0^ 
+  9^(5)  cosBjX  +  36n(5)  Sinf.3x 


c.,x  -a,x 

■*  97n(5>e  +  98„<?>S 


8 

^  Sjn^VjW 

i=i 


026) 
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and  similarly: 


where: 


and: 


_a  p 

gln(§)  =e  5  (K1  cos35?  -  K2  sin05§) 

-a.?  _  _ 

g2n(?)  =  e  (K]  sing5§  +  K2  cos65?; 

&C?  _  _ 

93n(?)  =  e  (K 1  cosP5?  +  K2  sin05§) 

a-c?  _  _ 

9^(5)  -  e  (Kj  sin05?  -  K2  cosp5?) 

n-1 

95n^)  ”  (-1)  2  K203COS¥ 
n-1 

g6n(l)  =  H)  2  K3p"sin03§ 

1  r,  "al? 

97n(§)  =  2alKle 

tc\  1  a 

98n(?)  =  2alKle 


-a..?  _  _ 

glm(?)  =  e  (_K3  sing5?  ~  K4  cosP5?) 

-ct5?  _ 

g2mte)  =  e  (^3  cosg^  -  K4  sinp5§) 

a5?  _  _ 

g3m(?)  =  e  (_K3  s5n|35?  +  K4  COs35?) 

g4m(=)  =  e  5  (K3  cos^5^  +  K4  sin^ 


g5m(?)  =  -H)2  K2P3  sinrf^ 
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m 


g6m{?)  =  (_ir  K2^3COS¥ 


97rrP  2alKle 


ap 


z-\  _  *  ni  1/  1 

98m(?,“2al  kle 

K.  =l(a  K,  +  b  K  ) 

1  2  n  3  n  4 

K_  =  i(a  K,  -  b  KJ 

2  2  n  4  n  3 

K.  =  i(a  K_  +  b  KJ 

3  2  m  3  m  4 

K.  =.i  (a  K.  -b  KJ 

4  2  m  4  m3 


In  addition  to  (124)  the  following  conditions  must  be  satisfied: 

j* 

•— r L(G)  - 0  i  =  2,3,--.  (127) 

dx' 

where  L  is  the  differential  operator  defined  earlier.  Equations  (127)  lead  to  the 
expressions: 


9  C 

dVG(0)  _  ^2 


dx 


_2 

C1 


d]lG(0)..  i  ,2  c  , 

,11  ~„3V^2  T3J 

dx  C  ^ 

d'3G(0)  _  1  c ■  r  r2„  3 

-7T3~-^<2ClC2C3'C^4-(-2) 

dx  C  j 


The  solution  (121)  may  be  written  in  the  form: 

8 


o 

7 

xz 


-  V  [A.+  A.(x)3s(x) 

- — r  i  f  i 

i=i 


028) 


029) 
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where: 


7  x 

'|W=I  /  Fk«»jk«« 

k~2  o 


it  can  be  shown  that  the  transverse  shear  force,  V  ,  the  stress  resultant.  N  ,  and  the 

xu  xu 

stress  couple,  M  ,  can  be  expressed  as: 


V  =  V  [A.  +  A ;(x)]$ .(*) 


N  = 
xu 


=  S  t  V  A  (x)i  (x)  +  A“  jAq  +  B“, A10  -  N;u 


where  the  functions  ty.(x),  pL.(x),  and  p.(x)  have  been  defined  as  part  of  equations  (96). 

The  average  slope  cf  the  lower  laminate  at  x  -  0  and  the  average  slope  of  the  upper 
laminate  at  x  =  L  are  given  by: 


w'x  i  i  (P4AV  "  R22)cpj^°)  +  (P6AV  "  Z22?*j'HAj 

i=1 

+  A!l  +  Cj"  1^22  *  P4AV)(Z22ri2  "  Z12V+ 

(P6AV  "  ^22)(R22T12  "  R*2n4)l  "  2  a*1  "  Ab2*  Sp,> 
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«'x  ■  S  Kw  (PMV  -  Z22)*i<L)|lAi  +  Ai<L>l 

1  =  1 

+  a]0l  +  A'n  +  |(»  /2  ^  w)(z*2n2  -  z*2n4) 

+  *P6AV  "  Z22)^R22T12  "  R1 2V I  +  2  fa(1  +  Ab2*  ezp,x  *132) 

where  Pb  ,,  PA.W,  Pb  ,,  and  P^AW  are  as  given  inequation  (91).  Instead  of  equation 
4AV  4AV  6AV  6AV  0  n  t 

(96)  one  has  for  the  non-linear  problem: 

[H]{A*}  =  {C*}  +  (Ap)  (133) 

which  differs  from  the  linear  case  only  in  the  addition  of  the  column  matrix  {A^} .  This 
matrix  contains  all  the  plastic  strain  terms.  The  non-zero  coefficients  of  {A^}  may  be 
obtained  by  using  equations  (129),  (130),  (131),  and  (132)  in  place  of  the  corresponding 
elastic  relations,  in  the  formulation  of  the  boundary  conditions.  They  are  given  by: 


V2)  =  ■  S  Ai(L)cDi{L) 

Pi  ‘ 


A  (4)  =  -  S  A.(LH.(D 

P  ^  I  I 


A  (6)  =  -  y  A.(L)U  (L) 

P  I  I 


V9)  =  - 1  Ai(L)pi(L) 
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V,0>  =  I  Ai(L)  I  (if22  -  P4AvK<L>  *  (Z22  -  pU6AV>*i'<l>l 


i=i 


1  I/2T 


|(R22  '  P4AV)*Z22t,2  "  Z*2T,4)  *  *Z22  "  P6Av'*R22T,2  "R12,,4)l 


yn>=- 


+  7T  t  (1  +  AbJ 
2  a  2 


C]  !(R22  "  P4AV)(Z22ri2  “  Z*1 2n4)  “  ^Z22  “  P6AV^R22Ti2  "R12T14)I 


zp,x 


2  ta(l  ~^2] 


zp,x 


In  order  to  determine  the  coefficients,  A?,  in  equations  (133),  the  coefficients,  must 

be  known.  This  in  turn  requires  a  knowledge  of  the  plastic  strain  distributions  in  the 

bond-layer.  The  solution  to  this  problem  may  be  obtained  by  successive  approximations, 

starting  with  the  assumption  that  the  plastic  strains,  and  hence  the  coefficients  of  (A  3  / 

P 

are  zero.  Initial  values  for  the  plastic  strains  are  therefore  obtained  based  on  elastic 

adhesive  stresses.  The  piastic  strains  are  computed  at  a  number  of  stations  along  the 

length  of  the  joint,  so  that  a  numerical  integration  procedure  may  be  used  to  evaluate  the 

integrals  necessary  to  determine  the  coefficients  A.(x),  defined  in  equation  (129).  After 

calculating  the  coefficients  of  {A  },  new  values  for  A*  may  be  obtained  by  solution  of 

P  I 

equations  (133).  This  leads  to  a  new  set  of  adhesive  stresses  and  plastic  strains.  The 
above  process  is  repeated  until  the  desired  convergence  is  obtained. 


Although  the  procedure  described  above  appears  sound,  a  number  of  computational 
problems  dealing  with  accuracy  and  convergence  are  anticipated.  Since  the  development 
of  a  non-linear  analysis  program  of  this  complexity  is  outside  the  scope  of  the  present 
contract,  it  was  decided  to  develop  an  approximate  but  far  less  complex  approach  as  an 
interim  method  of  analysis. 
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4.5.4  Plastic  Zone  Approach 

An  analysis  procedure  and  associated  computer  program  was  developed  for  bonded  joints 
with  non-linear  adhesive  stress-strain  behavior.  The  analysis  uses  a  "plastic  zone" 
approach  together  with  a  two-stage,  elasto-plastic,  effective  stress -stra in  curve.  The 
assumption  that  plastic  zones  develop  at  the  ends  of  the  over -lap  when  the  adhesive 
stresses  reach  a  certain  (maximum)  value,  and  then  spread  inwardly  toward  the  center  of 
the  overlap  when  the  load  on  the  joint  increases,  was  reported  by  Douglas  Aircraft 
Company,  Inc.  (referenced).  The  Douglas  analysis,  however,  ignores  the  presence  of 
normal  stresses  in  the  bond-layer  and  assumes  that  the  shear  stresses  are  constant  inside  the 
plastic  zones.  In  practice  neither  the  shear  stress  nor  the  affective  stress  will  be  constant 
for  any  appreciable  distance  along  the  bond-layer.  Nevertheless,  it  was  felt  that  inte¬ 
grating  the  "plastic  zone"  concept  with  the  linear  bonded  joint  analysis  (BONJO  I)  would 
provide  a  good  approximation  for  the  analysis  of  bonded  joints  ir,  which  the  adhesive  was 
stressed  above  its  yield  strength. 

The  effective  stress  is  defined  accoring  to  the  Von  Mises  condition  as: 

_ 

a  =  \  o°  +3 a° 

V  7  xz 

which  is  identical  to  the  first  of  equations  (113).  The  effective  stress  inside  the  plastic 
zones  is  assumed  constant  and  equal  to  the  maximum  stress  obtained  from  a  unidirectional 
stress-strain  curve.  The  stress -strain  behavior  of  the  adhesive  between  the  plastic  zones  is 
assumed  to  be  linear  elastic;  hence,  the  slope  of  the  effective  stress-strain  curve  in  this 
region  will  be  equal  to  3G .  As  illustrated  in  Figure  17,  the  two-stage  stress-strain  curve 
used  in  this  procedure  is  obtained  by  extending  the  elastic  region  until  the  ultimate  stress 
is  reached,  at  which  point  the  strain  increases  without  further  increase  in  stress.  Deter¬ 
mining  the  length  of  the  plastic  zones  is  a  quickly  convergent  iterative  process.  Initially 
an  elastic  solution  is  performed  and  adhesive  shear  and  normal  stresses  are  calculated  at 
76  locations  as  indicated  in  Figure  18.  The  initial  lengths  of  the  plastic  zones  are 
selected  so  that  the  effective  stresses  are  below  the  maximum  stress  level  at  all  locations 
situated  between  the  two  zones.  After  the  plastic  zone  lengths  are  obtained,  the  shear 
and  normal  stresses  within  those  zones  are  reduced  proportionally  so  that  the  resulting 


81 


FIGURE  17  EFFECTIVE  STRESS -STRAIN  CURVE  EPON  96OI 
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FIGURE  18  EFFECTIVE  STRESSES  -  PLASTIC  ZONE  APPROACH 


effective  stress  becomes  equal  to  its  maximum  value,  or: 


a 

o  max  o 

a  = - a 

xzp  -  xz 

r  a 


a 

max  o 

- a 

—  z 
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034) 


In  the  next  step,  an  elastic  analysis  is  performed  for  that  portion  ot  the  joint  which  is 
between  the  two  plastic  zones.  Although  the  analysis  itself  is  identical  to  that  discussed 
in  Sections  4.2  and  4.3,  new  boundary  conditions  must  be  established  for  this  problem, 
which  may  be  obtained  from  the  known  shear  and  normal  stresses  inside  the  plastic  zones. 
Denoting  the  length  of  the  plastic  zone  near  the  end  of  the  splice  plate  (x  =  0)  by  x^  and 
that  near  tne  center  of  the  splice  plate  by  X£,  the  following  boundary  conditions  are 
defined  for  the  center  portion  of  the  joint: 


035) 
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For  a  single  lap  joint,  the  condition 

M  (x)  M  .  (x)  =  M  1  36) 

xu  xL  x 

is  also  needed  in  order  to  obtain  a  solution.  The  total  moment,  M  ,  may  be  determined 
from  the  condition  that  the  slope  is  zero  at  the  ends  and  at  the  center  of  the  joint.  All 
integrations  in  equations  (135)  are  performed  numerically  in  the  program. 


The  effective  stresses  resulting  from  the  above  analysis  are  generally  as  shown  by  the 
dashed  line  in  Figure  18.  Since  the  effective  stresses  at  points  A  and  9  exceed  the  maxi¬ 
mum  stress  level,  the  plastic  zone  lengths  are  increased  and  the  entire  process  is  repeated 
until  the  difference  between  the  plastic  zone  lengths  of  successive  iterations  becomes 
negligible. 


Figure  19  shows  the  peak  adhesive  shear  stresses  and  the  final  plastic  zone  length  at  each 
side  of  the  joint  as  a  function  of  the  total  applied  load,  for  the  configuration  D  specimen. 
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FIGURE  19  PEAK  ADHESIVE  SHEAR  STRESSES  AND  PLASTIC  ZONE  LENGTHS 
NON-LINEAR  BONDED  JOINT  ANALYSIS 
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IV.  FINITE  ELEMENT  ANALYSIS 


1.0  GENERAL 

Finite  element  analyses  were  performed  as  part  of  this  contract  for  the  purpose  of 
evaluating  the  numerical  results  obtained  with  the  closed  form  analysis  program,  BONJO  I . 
These  results  were  presented  and  discussed  in  Chapter  HI,  Section  4.4.  In  addition, 
finite  element  cnafyses  were  made  for  some  of  the  step  joint  ana  mechanical  joint  speci¬ 
mens,  which  will  be  described  in  this  chapter. 

Lockheed  FAMAS  Program  ^97  was  uced  to  perform  the  numerical  calculations.  This  pro- 
giam  employs  the  direct-stiffness  displacement  method  to  perform  o  linear  structural 
analysis  for  deflecrioris  and  internal  loads  of  statical!/  loaded  structures.  Formulation  and 
decomposition  of  the  structural -stiffness  matrix  are  done  in  double  precision  .  Anisotropic 
h  iangular  constant  stress  elements  were  incorporated  in  this  program.  Several  options  ore 
available  to  input  the  material  properties  so  that  either  a  plane  stress  or  plane  strain 
analysis  can  be  made.  A  capability  to  determine  thermally  induced  stresses  or  strains  in 
anisotropic  structures  fas  recently  been  added  to  the  piogrom. 


2.0  STEP  LAP  JOINT  ANALYSES 

Analyses  were  mode  of  two  different  step  lap  joints  of  the  Configuration  'B'  small  scole 
specimens.  Detailed  dimensions  and  materia!  specifications  are  given  on  Drawing  No. 
7226~,302l3  in  Appendix  C  of  Volume  II.  The  first  joint  analyzed  utilizes  a  16-ply  boron 
Icminate  (-13  specimen)  consisting  of  8  plies  at  0°  and  8  plies  at  ±45°.  The  laminate  is 
bondeJ  to  an  aluminum  adherend  in  three  steps  of  0.50  inch  overlap  each.  The  adhesive 
thickness  was  taken  as  0.0042  inch,  which  was  the  approximate  overage  thickness  of  all 
test  specimens.  The  pertinent  port  of  the  finite  element  model  used  for  the  analysis  is  shown 
in  Figure  20.  A  finer  grid  is  used  near  the  ends  of  each  step  since  higher  gradiaits  exist 
in  these  areas.  Figures  21  and  22  show  the  adhesive  shear  and  normal  stresses,  respectively, 
along  the  length  of  the  join*.  Both  types  of  stresses  peal:  at  the  ends  of  each  step  and 
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FIGURE  20  FINITE  ELEMENT  MODEL-STEP  JOINT 


FIGURE  22  ADHESIVE  NORMAL  STRESS  -  STEP  JOINT 
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hove  their  maximum  volues  at  the  end  of  the  last  step  where  the  aluminum  to  boron  thick¬ 
ness  ratio  is  largest.  The  longitudinal  stresses  in  the  aluminum  adherend  along  three 
parallel  surfaces  are  presented  in  Figures  23.  These  surfaces  are  situated  so  that  each  one 
represents  the  mid-plane  of  a  row  of  elements  adjacent  to  the  bond-layer  of  one  of  the 
steps.  The  stresses  along  a  g  >en  surface  adjacent  to  the  bond-layer  build  up  rapidly 
toward  the  end  of  the  step  and  then  drop  sharply  at  the  beginning  of  the  next  step  when 
the  adherend  thickness  is  increased.  The  longitudinal  stresses  in  the  boron  laminate  are 
presented  in  a  similar  manner  in  Figure  24. 

A  modified  analysis  was  performed  to  investigate  the  effect  of  recognizing  the  tension  ties 
between  the  adherends  and  between  the  adhesive  and  adherends  at  the  ends  of  each  step 
The  resulting  adhesive  shear  stress  distribution  is  plotted  in  Figure  25.  For  comparison, 
the  adhesive  shear  stress  distribution  obtained  in  the  original  analysis  (without  tension 
ties)  is  shown  also.  It  can  be  seen  that  the  peak  shear  stresses  at  the  ends  of  the  steps  are 
reduced  significantly  as  a  result  of  recognizing  the  tension  ties  and  that  the  maximum 
adhesive  shear  stress  in  the  joint  was  cut  by  30  percent.  Because  of  the  fact  that  part  of 
the  applies  load  is  reacted  by  tensile  forces  between  the  adherends,  the  total  area  under 
the  shear  curve  no  longer  represents  the  applied  load  for  the  modified  analysis.  A  detail, 
showing  how  the  joint  was  modeled  at  the  ends  of  each  step,  is  presented  in  Figure  25. 

The  tensile  stresses  in  the  shaded  elements  range  from  2000  to  3000  psi,  while  the  shear 
stresses  in  these  elements  are  almost  negligible. 

The  second  joint  analyzed  was  identical  to  the  first  one  except  that  titanijm  was  substi¬ 
tuted  for  aluminum  for  one  of  the  adherends.  The  finite  element  model  shov/n  ir.  Figure  20 
therefore  remains  applicable.  Tension  ties  were  assumed  to  exist  at  the  ends  of  each  step. 
The  adhesive  shear  stress  distribution  for  this  case  is  shown  in  Figure  26.  The  corresponding 
stresses  for  the  case  with  aluminum  adherends  are  given  by  the  dashed  line,  for  comparison. 
The  bcron/titanium  joint  yields  considerably  lower  shear  stresses  at  the  end  of  the  last 
step  than  the  boron/aluminum  one,  but  generally  the  stresses  differ  only  slightly. 
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With  Tension  Ties 


FIGURE  2f.  ADHESIVE  SHEAR  STRESS  DISTRIBUTION  -  STEP  JOINT 
COMPARISON  OF  ALUMINUM  AND  TITANIUM  ADHEREND6 


3.0  MECHANICAL  JOINT  ANALYSES 


Two  separate  finite  element  models  were  constructed  in  order  to  perform  a  detailed 
analysis  of  one  of  the  Configuration  *E *  small  scale  specimens.  Dimsnsions  and  material 
specifications  are  given  on  Drawing  No.  7226-1 302IE  in  Appendix  C  of  Volume  II .  The 
-1A  specimen  assembly  consisting  of  an  8-ply  boron  laminate  (-21  specimen)  and  a 
titanium  splice  plate  was  selected.  Two  0.012  inch  thick  titanium  shims  were  inter- 
layered  with  the  boron  in  order  to  provide  su*ficient  bearing  for  the  fasteners.  The  first 
model  which  is  shown  in  Figure  27  is  for  the  purpose  of  performing  a  detailed  analysis 
through  the  thickness  of  the  joint.  Isotropic  triangular  plate  elements  are  used  for  the 
titanium  splice  plate/  titanium  shims,  adhesive  layers,  and  fasteners.  Anisotropic 
triangular  plate  elements  are  used  for  the  boron  laminae.  In  order  to  properly  account  for 
pin  bending  and  *o  determine  the  bolt  bearing  loads  on  the  individual  layers  of  the  joint, 
the  elements  of  the  fasteners  are  connected  to  those  of  the  joint  plates  with  springs  that 
are  permitted  to  take  compression  loads  only. 

In  order  to  study  the  stress  distribution  around  the  fastener  holes  a  second  finite  element 
model  was  constructed.  The  latter  model  consists  of  triangular  plate  elements  in  the 
plane  of  the  joint  and  represents  a  layer  of  titanium  or  boron  ,  Again,  compression  springs 
between  the  fastener  and  plate  elements  are  used  to  obtain  the  bearing  stresses  caused  by 
the  bolt  in  the  iayer  under  consideration.  This  second  model  has  been  shown  in  Figure  28. 

Average  bearing  stresses  in  the  two  fasteners  were  calculated  for  an  applied  compression 
load  of  1000^.  These  stresses  were  obtained  from  the  first  model  by  dividing  the  final 
forces  in  the  springs,  between  the  fasteners  and  joint  plates,  by  their  respective  areas  of 
contact.  An  iterative  analysis  was  used  to  assure  that  all  final  spring  forces  were  com¬ 
pression  by  successive  elimination  of  the  tension  springs.  The  average  bearing  stresses  thus 
obtained  are  shown  in  Figure  29.  Peripheral  bearing  stresses  on  the  fasteners  were  deter¬ 
mined  with  the  use  of  the  second  finite  element  model.  Figure  30  shows  those  transmitted 
by  the  upper  titanium  shim  as  a  function  of  the  angle  8. 
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FIGURE  27  FINITE  ELEMENT  MODEL  NO.  1  -  MECHANICAL  JOIN' 


TITANIUM  SHIM 


Loading  and  unloading  of  the  titanium  shims  has  been  studied.  Figure  31  shows  the  shear 
stresses  in  the  adhesive  layer  on  both  sides  of  the  upper  titanium  shim  resulting  from  an 
applied  compression  load  of  1000*.  Peak  shear  stresses,  of  course,  occur  near  the 
fasteners  where  tne  loads  are  introduced  into  the  shim.  The  position  of  the  fasteners  is 
indicated  in  the  figure. 


Net  section  tensile  stresses  for  the  lower  titanium  shim  and  the  0  boron  layer  in  the 
center  of  the  plate  are  shown  in  Figures  32  and  33,  respectively,  for  a  load  of  1000  lbs. 
These  stress'  /ere  determined  by  taking  the  average  stresses  acting  on  the  layers,  as 
obtained  from  n  first  P  ;te  element  modle,  end  applying  them  to  the  second  model . 
The  stress  concentraf’on  .-jetor  appears  to  be  significantly  higher  for  the  baron . 


FIGURE  31  SHEAR  TRANSFER  -  TITANIUM  SHIM 


FIGURE  33  NET  SECTION  STRESSES  -  O  BORON 


V.  PHOTOELASTIC  STRESS  ANALYSIS  -  BONDED  JOINTS 


1.0  GENERAL 

An  attempt  was  made  to  provide  an  experimental  verification  of  the  analytical  stress 
analysis  methodology  developed  under  the  program  (Chapter  III).  This  verification  was 
limited  to  the  linear  analysis  regime  because  of  its  relative  importance  in  serving  as  the 
foundation  of  the  subsequent  non-linear  analysis. 

Basically,  three  experimental  phases  were  accomplished  during  the  program.  The  first  of 
these  involved  the  fabrication  and  photoeiastic  test  of  a  Configuration  D  double  lap  joint 
having  aluminum  adherend  strips  and  glass  splice  plates.  This  phase  is  described  in 
Section  1.1.  During  the  second  experimental  phase,  Configuration  A  and  D  joints  were 
tested  to  determine  the  stress  distribution  at  the  surface  of  its  titanium  splice  plates. 

This  second  phase  is  described  in  Section  1 .2.  Finally,  Section  1 .3  describes  the  third 
experimental  phase  in  which  strain  gage  results  were  obtained  for  several  configuration 
A  and  D  specimens.  These  results  were  used  to  provide  an  independent  check  on  the 
ant  lytical  and/or  photoeiastic  results. 


1.1  PHOTOELASTIC  RESULTS  -  ISOTROPIC  ADHEREND 

This  phase  of  the  investigation  was  initiated  in  an  attempt  to  verify  an  existing  closed- 
form  solution  for  the  stress  distribution  in  a  linearly  elastic,  elastically  isotropic  joint 
adhesive.  In  order  to  achieve  this  goal,  a  photoeiastic  model  of  a  double  lap  joint  was 
constructed.  The  adhesive  layer  was  simulated  by  a  thin  layer  of  epoxy  photoeiastic 
material.  A  reflective  coating  was  painted  on  one  side  of  the  specimen,  between  the 
aluminum/epoxy  interface,  to  allow  determination  of  isochromatic  lines.  Although  the 
optical  sensitivity  of  glass  is  lov.  relative  to  that  of  the  adhesive  layer,  a  reflective 
coating  was  painted  between  the  epoxy/glass  interface  on  the  other  side  of  the  specimen 
to  determine  any  contribution  to  the  overall  fringe  pattern  made  by  the  gloss.  A  schemafi 
of  this  model  is  shown  in  Figure  34. 
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SECTION  A-A 


FIGURE  34  PHOTOELASTIC  MODEL 


The  model  was  constructed  to  allow  determination  of  the  stress  distribution  (0,0,0  ), 

x  y  xy 

in  the  plane  of  the  adhesive  layer.  Although  the  analysis  of  this  configuration  is  based  on 
the  assumption  that  these  stresses  are  zero  (for  relatively  "narrow"  joints),  it  was  antici¬ 
pated  that  at  least  some  sort  of  distribution  could  be  determined  photoelastically .  No  such 
distribution  could  be  obtained.  Two  models  of  this  type  were  constructed  and  subjected  to 
monotonicaily  increasing  sequ-vices  of  static  loads  until  the  glass  splice  plates  failed  in 
tension  (Figure  35). 


106 


GLASS  SPLICE 
PLATES 


APPLIED 

LOAD 


APPLIED  LOAD 


FIGURE  35  VIEW  SHOWING  FAILED  GLASS 

SPLICE  PLATES  ON  PHOTOELASTIC  MODEL 


In  both  cases,  the  epoxy  used  to  mode!  the  adhesive  layer  was  of  sufficient  sensitivity  so 

that  at  the  recorded  failure  load  level,  the  shear  stress,  a  ,  should  have  produced  at 

least  one  full  fringe  »n  the  central  portion  of  the  bonded  area.  The  fact  that  no 

appreciable  fringe  pattern  was  observed  in  this  area  when  the  model  was  viewed  normally, 

means  simply  that  the  stresses  a  ,  a  ,  a  were  of  much  lower  magnitudes  than  was  a 

x  y  xy  xz 

This  observation  is,  of  course,  in  direct  substantiation  of  the  assumption  made  in  ?he 

analysis  that  a  ,  a  ,  and  a  can  indeed  be  neglected. 

'  x  y  xy 


Although  the  anticipated  measurement  of  a  detailed  stress  distribution  was  not  accomplished, 

the  results  which  were  obtained  do  seem  to  achieve  the  desired  end  goal .  They  substantiate 

the  validity  of  the  assumption  that  the  adhesive  stresses  a  and  a  are  negligible  in  com- 

x  xy 

parison  with  a  and  a  . 

z  xz 
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1 .2  PHOTOSTRESS  RESULTS  -  COMPOSITE  ADHEREND 


During  this  phase  of  the  program,  an  attempt  was  made  to  obtain  a  correlation  between 
analytical  joint-stress  predictions  and  experimentally  determined  joint-stress.  The  experi¬ 
mental  information  was  obtained  by  bonding  pads  of  birefringent  material  to  the  splice 
plate  surfaces  of  several  joint  specimens. 

hree  such  specimens  were  considered:  a  1  .l"-wide  specimen,  a  3.0"-wide  specimen, 
and  one  having  a  width  of  9.0  inches.  All  three  types  of  specimens  had  0°/±45°  boron- 
apoxy  adherend  laminates  and  splice  plates  made  of  Ti-^I-4V  annealed  material .  The 
data  obtained  were  reduced,  using  the  shear  difference  method,  to  obtain  stress  and/or 
strain  components  at  several  desired  locations. 

1.2.1  1 .1  "-Wide  Specimen 

A  photograph  of  this  specimen  is  shown  below  in  Figure  36.  The  birefringent  pads  used  on 
this  specimen  were  cut  from  a  0.042"-thick  sheet  of  Photostress,  Inc.,  S— 1 6  material. 

The  titanium  splice  plate  was  0.039"  thick  and  the  base  adherend  was  cut  from  an  8-ply 
symmetrical  laminate. 

Two  models  of  the  type  shown  in  Figure  36  were  tested.  The  second  of  these  models  was 
tested  in  an  attempt  to  resolve  difference  between  experimental  and  analytical  results 
discovered  during  consideration  of  the  first  model. 

The  experimental  investigation  on  the  first  model  was  conducted  by  loading  the  specimen 
to  a  1000  Ibj.  load  level  and  using  a  reflective  polariscope  to  obtain  isochormatic  and 
isoclinic  photoelastic  data.  These  data  were  reduced  to  obtain  the  axial  stress  com¬ 
ponent,  a  ,  plotted  in  Figure  37.  The  1  KIP  load  applied  to  the  model  was  net  sufficient 
to  cause  appreciable  yielding  of  any  of  its  constituent  materials.  The  results  obtained 
were  therefore  compared  with  analytical  predictions  obtained  from  the  linear,  closed-form 
analysis  developed  under  this  program. 


PHOTOELASTIC  COATING  OVER 
TI  SPLICE  PLATE 


FIGURE  36  1.1" -WIDE  CONFIGURATION  "D”  SPECIMEN  WITH 

BIREFRINGENT  PAD  BONDED  TO  SPLICE  PLATE 

The  analytical  results  obtained  from  the  closed-form  analysis  do  not  agree  with  photostress 
results  within  acceptable  tolerance  (Figure  37).  Although  the  shapes  of  the  "c losed-form" 
and  "photostress"  curves  are  similar,  the  difference  in  stress  magnitudes  is  thought  to  be 
too  great  to  ignore.  A  complementing  finite  element  analysis  was  done  in  an  attempt  to 
help  resolve  this  difference.  The  result  of  this  analysis  is  also  shown  in  Figure  37.  This 
latter  (finite  element)  analysis  tends  to  lend  more  credence  to  the  "closed -form" 
analytical  results  than  to  the  experimental  results. 

The  experimental  curve  in  Figure  37  was  corrected  for  tensile  reinforcing  effect  arising 
from  the  fact  that  the  coating  material  itself  carries  some  load.  The  derivation  of  this 
correction  factor  is  based  on  the  assumption  that  the  coating  is  in  a  constant  uniaxial 
state  of  stress.  This  assumption  is  obviously  incorrect  since  the  splice  plate  undergoes 
bending  as  well  as  longitudinal  tension.  A  bending  correction  should  therefore  be  applied 
to  the  photostress  results.  Since  the  degree  of  splice  plate  bending  cannot  readily  be 
assessed,  no  such  bending  correction  can  be  made. 
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Distance  from  End  of  Splice  Plate,  x 
Length  of  Splice  Plate  L 


FIGURE  37  -  COMPARISON  OF  ANALYTICAL  WITH  EXPERIMENTAL  SPLICE  PLATE  SURFACE 
STRESS  -  FIRST  1.1"  WIDE  CONFIGURATION  "D"  SPECIMEN 


An  attempt-  was  made  to  further  resolve  the  differences  between  the  three  curves  shown  in 
Figure  37.  This  was  done  by  placing  strain  gages  at  selected  locations  along  the  center- 
line  of  the  surface  of  the  splice  plate  of  the  second  1  .  l"-wide  joint.  The  second  joint 
specimen  was  a  1 6— ply  thick  laminate.  After  strain  gage  results  were  obtained,  the  gages 
were  removed  and  photostress  material  bonded  to  the  splice  plate  surfaces.  A  photostress 
investigation  was  then  performed  on  the  same  model.  These  new  strain  gage  and  photo¬ 
stress  results  were  compared  not  only  with  each  other,  but  with  closed  form  analytical 
results  as  well .  The  results  of  this  comparison  are  presented  on  Figure  38.  It  can  be  seen 
that  these  three  types  of  results  correlate  much  better  than  comparable  results  obtained  from 
the  first  specimen . 

1.2.2  3.0“-Wide  Specimen 

Figure  39  shows  a  sketch  of  the  3.0"-wide  specimen  discussed  in  this  section.  The  speci¬ 
men  was  of  the  3.0"-wide,  Configuration  "A"  variety.  The  0°/±45°  base  adherend 
material  was  fabricated  of  Narmco  5505  boron/epoxy  and  wjs  an  8-ply  thick  symmetrical 
laminate.  The  splice  plate  was  cut  from  a  sheet  of  Ti-6AI-*'V  annealed  material  and  had 
a  0.0425"-thick  birefringent  pad  bonded  to  a  portion  of  jurface  as  shown  in  Figure  39. 

Photostress  data  obtained  at  each  of  the  grid  nodes  were  used  to  obtain  normalized  plots 
of  splice  plate  surface  stress.  Figures  40  and  41  show  the  longitudinal  and  transverse 
stress  components,  respectively.  The  longitudinal  component  displayed  in  Figure  40  was 
evaluated  along  line  ft)  (Figure  39) .  The  transverse  component  of  Figure  41  was  evaluated 
along  line  AB. 

The  obvious  anticipated  symmetry  about  lines  AB  and  CD  is  not  particularly  evident  from 
the  experimental  results  presented  in  Figures  40  and  41 .  This  would  seem  to  indicate  that 
the  data  on  which  these  results  are  based  are  somewhat  questionable.  The  explanation  for 
this  apparent  discrepancy  may  well  be  the  bending  undergone  by  the  splice  plate  as  a 
result  of  the  asymmetry  of  the  joint  about  its  midplone.  In  order  to  properly  interpret  the 
photostress  data,  a  correction  factor  must  be  used  to  account  for  the  degree  of  bending 
caused  by  this  asy  tmetry.  Such  a  factor  cannot  readily  be  obtained  for  the  situation 
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FIGURE  40  SPLICE  PLATE  SURFACE  LONGITUDINAL  STRESS  COMPONENT - 
3.0"  WIDE  CONFIGURATION  "A"  SPECB€N 
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herein  discussed  and,  therefore,  no  such  correction  was  made  to  the  results  presented  on 
Figures  40  <md  41 . 


1.2.3  9.0“ -Wide  Specimen 

Figure  42  shows  a  sketch  of  the  9.0"-wide  specimen  discussed  in  this  section.  The  speci¬ 
men  was  of  the  9.0"-wide,  Configuration  “A"  variety,  having  a  base  adherend  fabricated 
of  Narmco  5505  boron/epoxy.  The  laminate  was  of  an  8-ply  thick,  {0°/±A5°/0°)2  con¬ 
figuration.  Three  pads  of  b-refringent  material  were  bonded  to  the  splice  plate  at  the 
three  locations  shown  in  Figure  42. 

Photostress  data  obtained  at  each  of  the  grid  nodes  were  used  to  obtain  normalized  plots  of 
splice  plate  surface  stress.  Figures  43  and  44  show  the  longitudinal  and  transverse  stress 
components,  respectively,  evaluated  along  selected  lines  on  each  of  the  three  photostress 
pads. 

In  similarity  with  results  presented  previously  for  the  3.0" -wide  specimen,  the  symmetry 
about  the  specimen  centerlines  is  not  particularly  evident  from  the  plots  of  Figures  43  and 
44.  It  is  thought  significant,  however,  that  the  shapes  of  the  curves  displayed  on  Figures 
40  and  41  as  compared  with  those  displayed  on  Figures  43  and  44,  respectively,  are  similar. 


1.3  STRAIN  GAGE  RESULTS 

Strain  gage  results  were  obtained  on  ail  three  specimen  widths.  These  results  were 
obtained  by  placing  strain  gages  on  the  surface  of  the  titanium  splice  plate  to  obtain  data 
along  the  same  lines  at  which  photostress  data  were  obtained.  Each  such  cluster  of  gages 
was  arranged  as  shown  in  Figure  45.  Gages  1  through  5  were  oriented  long itudi nelly  and 
gages  6 and  7,  transversely. 

The  results  obtained  from  the  second  1  .l"-wide  configuration  "D"  specimen  are  plotted  on 
Figure  38.  These  results  are  compared  with  corresponding  photostress  and  analytical 
results.  The  maximum  joint  load  of  3000  lb^  was  chosen  to  avoid  yielding  of  either  splice 
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FIGURE  42  9.0" -WIDE  CONFIGURATION  "A"  SPECIMEN 

USED  TO  OBTAIN  PHOTOSTRESS  DATA 
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FIGURE  43  SPLICE  PLATE  SURFACE  LONGITUDINAL  STRESS 

COMPONENT -9*0"  WIDE  CONFIGURATION  "A”  SPECIMEN 
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FIGURE  45  ARRANGEMENT  OF  STRAIN  GAGES  j.N  AREA  OF  PHOTOSTRESS  PAD 


plate  or  adhesive.  Strain  readings  of  all  seven  gages  were  recorded  at  300  Ib^.  load  incre¬ 
ment  up  to  the  maximum  (3000  Ib^)  load  level.  Strain-lood  plots  obtained  from  these 
readings  were  judged  to  be  sufficiently  linear  to  satisfy  the  purpose  of  the  investigation. 

Gage  No.  6  exhibited  the  most  pronounced  deviation  from  linearity.  The  strain -load  plot 
for  this  gage  is  shown  in  Figure  46.  Figure  38  shows  that  the  c’osed  form  linear  analytical 
predictions  compare  favorably  with  strain  gage  results  along  the  longitudinal  centerline. 

Strain  gage  results  obtained  from  the  3.0"  wide  joint  specimen  are  displayed  in  Figures  47 
and  48.  Comparable  results  obtained  from  the  9.0"-wide  specimen  are  shown  in  Figures 
49  and  50 . 
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LONGITUDINAL  DISTANCE  FROM  SPLICE  PLATE  EDffl 

FIGURE  47  SPLICE  PLATE  SURFACE  LONGITUDINAL  STRAIN  COMPONENT  - 
3.0"  WIDE  CONFIGURATION  "A"  SPECIMEN 


TRANSVERSE  DISTANCE  FROM  CENTERLINE  OF  BIREFRINGENT  PAD 


FIGURE  48  SPLICE  PLATE  SURFACE  TRANSVERSE  STRAIN  COMPONENT  - 
3.0"  WIDE  CONFIGURATION  "A"  SPECIMEN 
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LONGITUDINAL  DISTANCE  FROM  SPLICE  PLATE  EDGE 


FIGURE  1+9  SPLICE  PLATE  SURFACE  LONGITUDINAL  STRAIN 

COMPONENT  -  9.C"  WIDE  CONFIGURATION  "A"  SPECIMEN 


TRANSVERSE  DISTANCE  FROM  CENTERLINE  OF  BIREFRINCENT  PAD 


FIGURE  50  SPLICE  PLATE  SURFACE  TRANSVERSE  STRAIN  COMPONENT  - 
9.0"  WIDE  CONFIGURATION  "A"  SPECIMEN 
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